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ABSTRACT 


This thesis describes the development of a 
minicomputer-controlled steady-state NMDR spectrometer 
system and its subsequent use to study the = es quad- 
rupole spin-lattice relaxation in monocrystalline 
NaNO,. 

The technique is an extension of that used by 
Pound (1950) and Hughes and Reed (1971), in which the 
change in amplitude of one component of a quadrupole- 
split NMR spectrum is measured when another component 
is severely saturated. Since the NMR signals are weak, 
we developed an efficient technique for measuring their 
amplitude, in which the first derivative of the absorp- 
tion signal is utilized to "lock" the spectrometer 
frequency to the resonance center, while the second 
derivative signal is used as a measure of the NMR am- 
plitude. To take full advantage of this novel technique, 
a mMinicomputer was used to control, in real time, many 
Subsidiary experimental parameters, such as the magne- 
tic field strength. The minicomputer was also used to 
collect and process the data. This enabled the relaxa- 
tion data to be acquired on an around-the-clock basis. 
To obtain meaningful results from such a system, great 


attention was paid to the elimination of systematic 


errors; these are discussed in detail in this thesis. 
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Extensive tests have shown that the system is capable 
of measuring NMR amplitudes LOraneaccuracysOL, about 
aU ad Ae 

The equipment was used to measure the orienta- 
tion dependence of Wo/W, of 23Na in NaNO, at room 
temperature. The quantities Wy and Wo, which are the 
probabilities of the Am= +1 and +2 nuclear magnetic 
transitions respectively, can be expressed as a linear 
combination of the components of a fourth-rank tensor, 
the so-called M-tensor. 

Excellent agreement was obtained between the 
data and the theoretical form in which it was assumed 
that Wy and Ws are governed by the same M-tensor. 

This showed that the spin-lattice relaxation was 

caused by a multiphonon interaction. An upper limit of 
eae ae was placed on the relaxation due to diffus- 
ing charges. The following M-component ratios were 


found: 03902542 000355 M 


3333/4111 = 
Ee On 1001 0003 LT eCone 


My 373/™@1311 = 
0.7091+ 0.0083, M 


bigot 
parison of these ratios with the predictions of a 

point-charge model favours the harmonic rather than 
the anharmonic Raman process. Comparison with the 
previous room temperature study made by Reed (1970) 
revealed significant discrepancies. [In particular, 


we found M123 to be zero, in agreement with crystal 


symmetry considerations but in disagreement with 
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Reed's results. We present possible explanations for 
the discrepancies. 

The time-dependent behaviour of the spin-system 
was experimentally studied in order to derive the value 
Of Wy directly. From these results we found that 
e*0*m Wels a0n0S00<s106 G4) Ne THiSevaiic 
is compared with that obtained by Niemel3 (1967) at 
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CHAPTER. I 
INTRODUCTION 


Nuclei with spin quantum number I > 1 can 
interact with their surroundings in two ways; a mag- 
netic interaction involving the nuclear magnetic 
dipole moment and the magnetic field at the nuclear 
Site, and an electrostatic interaction involving 
the nuclear electric quadrupole moment and the elec- 
tric field gradient at the nuclear site. The time- 
averaged magnetic and electrostatic interactions 
determine the structure of the nuclear magnetic 
resonance (NMR) spectrum. The so-called high field 
case is where the sample is situated in a large mag- 
netic field such that the Zeeman interaction is large 
compared with the static electric quadrupole inter- 
ScGtoOn. emormmcquivalent nuclei with spin number I = 
3/2 situated in a solid at sites of lower than cubic 
symmetry, one obtains a three-line spectrum 
provided the sample is in the form of a single crystal. 

The so-called spin-lattice interaction caused 
by time dependent magnetic fields and electric field 
gradients, results in a coupling between the nuclei 
and their surroundings (the lattice). The magnetic 
interaction can change the nuclear magnetic quantum 


number m by +1. In contrast, the electric quadrupole 
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interaction can change m=eby 21 or +2. The correspond- 
ing transition probabilities for the case LD = 3/72 are 
denoted by Wy and W, respectively. 

Pound (1950) performed a steady-state double 
resonance experiment on ONG Nuciewe (fom which gi 39/2) 
in sodium nitrate. By measuring the change in ampli- 
tude of one resonance when one of the others is 
Severely saturated by a large rf field, he was, able 
to demonstrate that the spin-lattice relaxation 
mechanism was quadrupolar rather than necrecics Later, 
Abragam (1961) pointed out that it was possible to 
obtain values of Wo/W, from a double resonance experi- 
ment of the type used by Pound. 

Using an extension of Pound's method, Reed (1970), 
in this laboratory, made a detailed Study of the orien- 
tation dependence of Wo/W, is@he 2 na ine sOdiumenitrate 
at room temperature, 

Pietila (1968) and Snyder and Hughes (1971) 
showed that Wy and W., are governed by a fourth-rank 
tensor, the so-called M-tensor. The sodium nuclei in 
sodium nitrate are situated at sites possessing 3 
Symmetry, and it can be shown that there are in this 
case only five independent M-tensor components. By 
fitting the orientation dependence of Wo/Wy to the 
theoretical form, Reed obtained ratios of the M-tensor 


components. His work indicated that the same M-tensor 
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governs Wy and Wo. This implies (Snyder and Hughes, 
1971) that the spin-lattice interaction is a Raman- 
type multiphonon interaction (Van Kranendonk, 1954; 
Van Kranendonk and Walker, 1967, 1968). Reed found 
generally good agreement with the theoretical form, 


except that the M component was non-zero, in 


L123 
disagreement with crystal symmetry considerations. 

In the meantime, Niemela (1967) measured the 
orientation dependence of Wy and Wo for the same 
nuclear spin system at 77 K. Niemela obtained his 
values of Wy and Wo by studying the approach to 
equilibrium of the spin system. His measurement 
technigue was different from Reed's in that the 
resonances were observed using a pulsed NMR technique 
(Hahn, 1950), rather than a steady-state method. 

Niemela unfortunately misinterpreted his data. However, 
the M-tensor values were later derived from Niemela's 
results by Hughes et al (1970) and were found to be 
rather similar to the room temperature values obtained 
by Reed (1970). Niemela's data did not allow the value 
of Mi 123 to be determined in his crystal. 

The work described in this thesis was undertaken 
in order to investigate further why the M1123 component 


was non-zero in sodium nitrate. A steady-state tech- 


nique was chosen in preference to a pulsed technique 
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Since the aoa resonances are rather closely spaced 

and thus could not be ee lecnet, excited by the 
pulses. However, to obtain an accuracy of 11% in 
Wo/Wy, it is necessary to measure resonance amplitudes 
to an accuracy of 70.1% in the steady-state double 
resonance technique. To obtain such accuracy obviously 
requires a sophisticated detection system in which 
particular attention is paid to linearity for example. 
Moreover, the signal-to-noise ratio is such that several 
days of data collection are required at anit crystal 
Orientation. It was thus felt that a computer con- 
trolled data acquisition system operating twenty-four 
hours per day was required. Reed's method consisted 

of recording the complete resonances in first deriva- 
tive form, and using the peak-to-peak amplitude as a 
measure of NMR signal amplitude. However, in this 
method, only a small fraction of the data collected 

is used. We therefore used a different method for 
measuring the NMR signal amplitude. In our technique 
which we believe to be novel, the first derivative 
Signal is used to lock the spectrometer to the resonance 
center, and the second derivative signal is used as a 
measure of the NMR signal amplitude. However, to take 
full advantage of this more efficient technique, 
tighter control has to be exercised over several 


experimental parameters. It was decided that the only 
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practical way of controlling such a complex system was 
by means of a real-time minicomputer. 

Considerable time was spent in designing and 
constructing the system. After the system was opera- 
tional, it was discovered that the data obtained 
suffered from a number of systematic errors. The 
elimination of these errors proved to be difficult, 
although illuminating. However, the final configura- 
tion seems to be capable of measuring resonance ampli- 
mudeSectO wan accuracy, Of about £0.12; 

A detailed study of the orientation dependence 
of W/W, of 2 NA in sodium nitrate at room temperature 
was then performed. Significant discrepancies were 
found between our results and those of Reed, and we 
present possible explanations for these. 

The steady-state double resonance technique does 


not give the magnitude of Wo Ol Woe but merely "gives 


2 
eheir ratio. sit is useful, and especially so for 
comparing measurements at different temperatures, to 
measure these quantities separately. The flexibility 
afforded by the minicomputer enabled the approach to 
equilibrium of the spin system to be experimentally 


studied. From this so-called transient experiment, 


the value of Wy was obtained. 
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The theory is presented in Chapter II and is 
applied to the specific case of e°NG in sodium nitrate 
rns Chapter itl. The method of obtaining W/W, 
and Wy from the experimental data is described in 
Chapter IV. The experimental details are presented 


in Chapter V, the results obtained are given in Chapter 


VI and are discussed in the following chapter. 
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CHAPTER OTT 
THEORY 


2.1 Basic theory of nuclear magnetic resonance 


We shall consider a simple description of the 
phenomenon of nuclear magnetic resonance. If we have 
a nucleus with a magnetic moment yu and colinear angular 
momentumetieimmersed in a magnetic field Hoe then the 


Hamiltonian of the system is 
tH nl be ean 1 (Zsa) 


where Ho Tsealongs them2 axis’ I, vse Cher zecompenents of 


I, and the gyromagnetic ratio y is defined by 


i a} / Flies ©. (222) 


This is the well known Zeeman interaction, and the eigen- 


values of H are given by 


En = ~yhH om (2) 


where the magnetic quantum number m can take the values 
py) vo. eee Le maximum component Olle rs) equadmrc 
Ih where the so-called nuclear spin quantum number I 
has either integral or half integral values. 

To detect such a system of energy levels one 
needs a resonant interaction that can cause transitions 


between them. In order to cause transitions between the 
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m and m' energy levels, a radiation field of angular 


frequency w must satisfy the condition 
Aw = ae EWM 2 (224) 


An alternating magnetic field perpendicular to Ho will 
cause magnetic transitions between energy levels rela- 
ted by m'=m+l. We therefore obtain the magnetic 


resonance condition 
w = yH : (255) 


The frequency w is also known as the Larmor frequency 
Since it coincides with the classical precession fre- 
quency of a magnetic moment in a steady magnetic field. 
For fields of order 10° gauss, the resonant frequency 
for nuclei occurs within the radio frequency (rf) 
Bange Ocleto SOsMHz. 

According to first order perturbation theory, the 
probability of these magnetic transitions is given by 
(Bloembergen, Purcell and Pound, 1948; hereafter referred 


to as BPP) 


2 
ede ty es Chore te (1/4) y7HSg (w) EPS) (GN epi ely ie (26) 


Here, Hy is the amplitude of the circularly polarized rf 
Magnetic field and g(w) is a normalized lineshape 
function which takes into account the finite width of 


the magnetic resonance signal. For an ensemble of 
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nuclei, in a solid for example, transitions with pro- 
bability given by (2.6) will cause the energy levels 

to become equally populated. The spin system is then 
said to be completely saturated, and its temperature 
gouinbinitesesin prackweenpetherspin system is embedded 
iieamelattice swith which 26 interacts. However, essen-— 
tially complete saturation can still be obtained provided 
this so-called spin-lattice interaction is weak compared 
with the interaction with the Hy field. In the absence 
of a Hy field, the spin system will come into equilibrium 
with the lattice, and the ratio of the populations of 
adjacent energy levels is given by the Boltzmann factor 
exp|AE/kT,|, where AE is He thveetgand aEipis: the lattice 
temperature. . Correspondingly, the ratio,of the -upward 
and downward spin-lattice transition probabilities is 
equal to the same poke naan Pact mar or nucle ata ttrcom 
temperature, AE is typically much less than kT, and the 
exponent can be replaced by 1+ (AE/kT,). 

If we now consider the case of I=1/2, eg protons, 
then the lower of the two energy levels will be slightly 
more populated than the upper. Thus when Hy is applied 
and the resonance condition is satisfied, there is a 
nets absorption ofvenergysifrom the iradiation field. tMThe 
magnetic resonance condition can therefore be detected 
by monitoring the radiation field. The field Hy will 


tend to equalize the energy level populations, whereas 
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the spin-lattice interaction will tend to maintain the 
spin system in equilibrium with the lattice. From the 
principle of detailed balancing, we find that the 

difference No between the population of the two energy 


levels satisfies the expression 


Ny == 2W(n,- N.) ~ 2NUP CP 
where from eq. (2.6), 
P= (1/4) y*H%g(w) (2.8) 


and Ny is the equilibrium value of Ne in the absence of 
H,- ined. eevee Wausethe Spin-latbice relaxation 

rate between the two energy levels. (Since the Boltzmann 
Pactore dat ters: so, lasttle from unity, there is no, need 

to distinguish here between the upward and downward 


transition rates.) In the absence of H the, solution 


LY 
Of) eg .46¢( 2.97.) «Us 


Nee noi - Cexp (-2Wt) ] C259) 


where C is a constant determined by the initial condi- 
tions fe:The timerconstane, of, the, decay. (onan is called 
the spin-lattice relaxation time T,- Ws ingmegs. (200) ) 
and (2.8), we find the steady state solution 


n 
SS eee ee 
beet 75 ; (210) 
b+ 3442) Hyg (w)T, 


an expression first derived by BPP. 
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The width and shape of the NMR absorption 
Signal is, in general, governed by several factors. 
Perhaps the most obvious is an inhomogeneity in Hi 
however, this is an experimental detail rather than a 
fundamental effect. The finite lifetime of the energy 
levels leads, according to the uncertainty principle, to 
a linewidth of order wee (expressed in frequency 
units). Another source of broadening is the magnetic 
dipolar coupling between nuclear spins. A typical 
nuclear magnetic dipole produces a field of order one 
gauss at a distance of three angstrom units. Thus in 
solids, the steady magnetic field will not be the same 
at each nucleus, but will vary by several gauss from one 
nucleus to another, thereby contributing to the*line- 
width. However, there is another contribution to the 
linewidth. One nucleus produces an alternating magnetic 
Prevasataithewstremor sanother., wilt hesnuclei ares identi-— 
cal, the alternating field is at the Larmor precession 
frequency of the second nucleus, and as a result may 
induce.a transition, This mutual interaction sreduces 
the lifetime of the energy levels of these nuclei, thus 
producing a "lifetime broadening", which turns out to be 
comparable to the other effect. This mutual exchange 
does not occur for unlike nuclei, whereas the first 


effect is still present. 
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The calculation of the NMR lineshape caused by 
Magnetic dipolar interactions is_a formidable task. 
However, aS was shown by Van Vleck (1948), the calcu- 
lation of the moments of the NMR line is much more 
tractable. The moments are usually expressed in fre- 


quency units, and are defined by 
M = | ( yg (w) d (aie 
se W- Wo) g(w)dw - 


where We is the resonance center frequency. In parti- 
cular, Van Vleck found that the second moment of the 
resonance of a spin species with gyromagnetic ratio Yr 


and spin I is given by 
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where Ys andes are the gyromagnetic ratio and spin of 
an unlike spin species, N is the number of I spins in 
eherunitrcel il, Box is the angle between Ho and the 

internuciearEvector Lak! and O50 Ss 


Hy and Si) s* ‘The summation j is over all spins I in 
7 ’ 


is the angle between 


the unit cell, the summation k is over all the I spins 
in the solid, and the summation & is over all the unlike 


spins S. 
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As indicated by BPP, another interaction between 
Nuclew andetnes Lactice to) possible, 31f) I > 1. In this 
case, the nucleus may possess an electric quadrupole 
moment which can interact with electric field gradients 
at the nuclear site. This interaction has three main 
effects; if the time-averaged electric field gradient 
is nonzero then additional structure in the NMR spectrum 
can result; these changes in the spectrum require modi- 
fications to the Van Vleck (1948) theory of second moment; 
Tietneselectric field gradient is fluctuating then an 
additional spin-lattice interaction is provided. These 
effects will be discussed in the following sections of 


this chapter. 


PeemeplLeccerr Cc Cuadeupole, moment. G7 the nucleus 


The Hamiltonian for the electrostatic interaction 


of the nucleus with electric fields is given by 
H = [p(x)v (x) av (22s) 


where p(r) is the charge density of the nucleus, V(r) 
is the electric potential produced by charges external 
to the nucleus, and the integral is taken over the 
nuclear volume. If we expand V(r) in a Taylor series 
about the nuclear center of mass, eq. (2.13) can be 


written in the following form 
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fe) ; 7 Poa 2 ca yk Poe, 
where 
2 es [xjo (dav 
and 


The first term in the expansion represents the electro- 
static energy of a point charge nucleus, and is thus 
independent of the nuclear shape or orientation. The 
second term is zero since the nuclear electric dipole 
moment ae is identically zero when measured relative to 
the nuclear center of mass. The third term involves 
the nuclear electric quadrupole moment tensor O44, and 
Phewerec trie field gradient tensar.,- The next term in 
the series involves the nuclear electric octopole 


Moment. However, this is again identically zero. . The 


following term involves the nuclear electric hexadecapole 


moment, which is also zero if I<2 (Abragam,1961; Cohen 
and Reif, 1957). Thus, in our work where we study the 
ona nucleus for which I is 3/2, we need only consider 
the electric quadrupole ianteraction. “Since the nucleus 
is in a state of definite angular momentum, i.e. the 
charge has cylindrical symmetry on account of its rapid 


spin, it can be shown (Abragam, 1961) that the orienta- 


tion dependent part of the nuclear electric quadrupole 
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interaction can be expressed in terms of a scalar quan- 
tity, the so-called electric quadrupole moment eQ which 


is defined by 
- we 2 
COPS mie eure TO (2515) 


where z* is. the axis of spin of the nucleus. 


2.3 Quadrupole split spectrum 


The Hamiltonian of a nucleus with quadrupole 
moment eQ which is situated in a magnetic field Hy and 
in an electric field gradient can be written (Abragam, 


1961) 
F e*qo 2 ih 2 2 
Hh = ~yhio- i +75 IT-1) [ate EET Lye Selsdlen, Pe ag A0G) 


where i 1, fil); Here the x, y, Z axes are .the prin- 
cipal axes of the electric field gradient tensor, 


arranged such that 


rvEe | tou (2027) 
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where vis is a°V/ax, ax THeSCLECL Ci Ceraeldscaradtenc 


B Ba 
eq and the asymmetry parameter n are then defined by 


the equations 


eq = V (2518) 
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(Since the electrostatic potential at the nuclear site 
satisfies Laplace equation y-V =.0, there are only 
two independent components Ver in the principal coor- 
dinate system. The quantities eq and n therefore 
completely specify the field gradient tensor.) 

If the site possesses spherical or cubic 
Symmetry, then by Laplace equation eg is zero. If 
the site is situated on an n-fold axis of symmetry with 
n 2 3, the asymmetry parameter is zero (Cohen and Reif, 
1957); eq is of course in general nonzero. For lower 
symmetry both eq and n are in general nonzero. 

The energy levels of the system can be readily 
found by a perturbation approach if the quadrupole 
interaction is either large (Das and Hahn, 1958) or 
small compared with the Zeeman interaction. In fields 
of order 16° Gauss, the electric quadrupole interaction 
of Se in solids falls into the latter category, the 
so-called’high field case. According to’ first’ order 
perturbation theory (Abragam, 1961), the energy levels 


are given by 


= il ae ~.S2 y 
En = hyom + (F)hvo- [3 COSs Us—-el+iSitim On cOSsse @ | 


-(m* - (3) 1(I+1)] (2.20) 


where the unperturbed frequency ve is equal to YH /2m 


and 
bo Ce, ieee epee (2-21) 
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The angles 6 and » are the polar and azimuthal angles 
of Hy relative to the principal axes of the electric 
field gradient. The second order perturbation term 


is in general quite complicated (Volkoff, 1953). Thus, 





we only present it for the special case of n = 0: 
vant 
(2) G8 Q 5 2 Lom (2 one. 
En = hla |n (3 cos §(1l-cos” 9) (8m~+1-4I (I4+1)) 
3 DIC 2 tn? 
= g (1-cos Oe 2m. 41-27 (T+1)) |. (252 2)) 


For the case ead Sy het O blows romeequs (2020) 
that the first order NMR spectrum consists of an unper- 
turbed center line corresponding to the m=5 <>-> 
transition, flanked by n pairs of satellite lines. Each 
member of a satellite pair is symmetrically disposed 
relative to the center line, and the maximum Splut eung 
is of order ge According to the second order perturba- 
tion theory, all the resonance meres are in general 
perturbed, the maximum shift being of order ee 
However, this shift is the same for both the upper and 


lower members of a satellite pair, so that their separa- 


tion is unaffected by second order quadrupole effects. 


2.4 Modifications to the dipolar broadening 


For the case of arbitrary spin and no quadrupole 
interaction, the resonance frequencies of the various 


magnetic transitions coincide. This degeneracy is 
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removed by the quadrupole shift in the resonance frequen- 
cies. Thus, the Van Vleck (1948) theory of dipolar 
broadening has to be modified in order to take into 
account the fact that not all nuclei now contribute to 
the broadening in the same way. The theory has been 
extended by Kambe and Ollom (1956) and Abragam (1961) 
for the center line resonance for half integer nuclear 
spin. 

Westiow detrne slike” nuclei to be nuciei of the 
Same species in equivalent sites having the same elec- 
meroms 1 eld gradien tl metThuss like nuctei have identical 


NRE opectra.) The contribution M. Oresuchenuciet to the 


2 
second moment of the center line is (Abragam, 1961) 
2 Zs eS 
Bie sear). (I+1) et Pea sao 
4 ny, rons ear epee \'t 
(1-3 BOP le Pe 
_) ———. (20922) 
fords 8 
Jr 3k 


Wheresthescummation 2S oOvervall “like” nuclei, and) tne 
remaining symbols have their usual meaning. 

If we define "semi-like" nuclei to be nuclei of 
the same species but ence at inequivalent sites in 
thesunit celle enone to first order, the center lines 
overlap and the satellite lines do not. The contribu- 
tion M, from "Semi-like" nuclei to the second moment of 


tne center line is (Abragam, 1961) 
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For unlike nuclei, i.e. those of a different 
Species, the second moment contribution is the same 
as that given by the Van Vleck formula (eqn 212) ))2 
Comparison of eqs.(2.23) and (2.24) with the Van Vleck 
formula shows that the center line is narrower than it 
is where there is no quadrupole Spliccing. 

According to Betsuyaku (1969), the second moment 
of the m «++ m-l resonance for the case of "like" spins 


is given by 





a Re al if 2 ome 2 2 - 2 
my = (ES + TTF 1 (I+1) L(I+1) (2a-1) + a“ +a 13] 
ie Gee cos.8r) 
-y,h a —— : CZ) 
j,k eae 
For example, it’ can be shown that for the case I = aye, 


the ratio of the second moment of a satellite to the 
second moment of the center line is 70.852 tfor thescase 


Of broadening by "like" sping. 
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Ales: Quadrupole spin-lattice relaxation 


We shail now consider the effect of a fluctuating 


electric field gradient at the nuclear site. Using the 


notation of Pietild (1968) and Snyder and Hines m og, 


the Hamiltonian of this interaction can be written in 


the form 
+2 aie 
Ho = ue Q (BY) (2.26) 
where 
Ons j 
BU = Vi ghe 
sot = 15 / Ga 
ul XZ pee ia 
is 1 Z . 
Bet (1/V6) (es ee ive) 


and 


ae = {312 SCI 
ot = tk (Y6/2) (1,1, + Peele) 
Q°* =k (6/2) 17 


SS Ay 


In a right handed coordinate system ey a cee 


which Ho lies along the z' axis (the x',y' axes being 


arbitrary) the matrix elements of Q" are 


I 


<mt+1|Q~*|m> 
ioe iy. ) 
+ (Y6/2)k (2m+1) { (Itm+1) (I¥m) } (2h 21.) 


<m+2 ise | m> 
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= (V¥6/2)k{(IFm) (I¢m-1) (T#m+1) (1+m+2) }27? (2.28) 


provided the quadrupole interaction is small compared 
with the Zeeman interaction. It can be shown by a 
perturbation approach that the relaxation probabilities 
for transitions which change the magnetic quantum number 


by one and two respectively are (Hughes, 1973) 


4 #3 
21d ) exp (“BE,) [<n" Vy - edie are 
ak db nn 
cee A eh Sg aoe Meet tuna: i. a. se 
y exp (-B8E_ ) 
nn' 
°d(E_,- E eee pe ) (22129) 
Ah nave eve ei eave cinsic 
25 cient Tes tact sah Seats 
Ww = = oy ge ee 
m>m+2 Fi y exp (-BE,) 
nn' 
SOA Eaie Bebe a cen Eis) (2.30) 
where 
Ay = (k/2) (2m+1) [ a (I+m+1) pve 
1/2 
Ao = (k/2) [ (I-m) (I-m+1) (I+m+1) (I+m+2) ] : 


Here n and n' are initial and final lattice states, En 
is an eigenvalue of the lattice and 8 is Cleans The 
m dependence of the multiplicative factors is removed 


by defining Wy and W, in the following way (Yosida ‘and 


Moriya, 1956) 


27. 
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W, (2m+1)* (I+m+1) (I=m) 








W = C2532) 
m>m+l1 21 (2T-1)7 
W. (I-m) (I-mt+1) (I+m+1) (I+m+2) 
W es ee (2832) 
m>m+2 21(2T-1)° 
Thus it follows that Wy and Wy can be written as 
s 2 
2.2 i; exp (-BE_ ) | <n NS ot [n> | 
w. = Ze Q nn' Y 
Tee T Tims) te eee ) : 
P n 
ee get ae ee ee) (27533) 
p BYE ms) a ae es: ! Z 
W me-07 ee aaa Hoh Seth Mine vi rig! LVyr yin? | 
= e Pe CS Ce oe Es 6 a A ee ee a eee 
a 4iii ) exp (-BE, ) 
nn' 
pA BR Et Ent2 ~ Em! Me EM 


Thesprobabilities Wy and Wo governethe relaxation, beha— 
viour of the various resonances in the quadrupole split 
spectrum. Their magnitudes are in general different 
because they involve different field gradient Components. 
However, they are expected to be of the same order of 


magnitude. 


2.6 Orientation dependence of the quadrupole relaxation 
a ee EE pO Lee relaxacion 


In the high field case, the external magnetic field 
Hy determines the axis of quantization. We have seen 
in section 2.3 that the NMR spectrum changes when we 


rotate Ho relative to the crystallographic axes. 
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Similarly, the quadrupole relaxation probabilities 
vary with the orientation of Hy relative to the 
crystal. 

To calculate the orientation dependence of Wy 
and W, wer transform, equations (2.33) and (2.34) Erom 
the primed coordinate system to an unprimed (right 
handed) coordinate system fixed with respect to the 
crystal. The polar and azimuthal angles of H, relative 
to the new coordinate system are denoted by 6 and 9. 
The transition probabilities Wi and W, can then be 


expressed as linear combinations of a fourth rank tensor 


Wye, 1957 Pietila, 1968) 


1 * 
) shee sol eal FS ete [¥gln> <n'|V [n> 6(E -E +E -E) 


L ieee NN my 
(Ret ey Se 
ab a's fi ) Exp (-6E, ) 
nn! : 
| (2.35) 
Wheres 6,0" 60 =) le 2.5 = Xeye Zz and q= 1) for Wy and 2 
for Wo- In general the Ligatg! are complex and thus we 
write 
Lygatg! = Migatg! + iN gatg? (2.36) 
from which it follows that 
Migate! = Mut grag (22.37) 
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Also, since Vagina Veo and from Laplace equation p Sb oen = NY, 


there are only fifteen independent M components 


apa'Bp! 
and ten independent Ni eitet components. 

If the lattice satisfies the principle of time 
reversal invariance (Goldberger and Watson, 1964) it 
can be shown (Hughes, 1973) that the quadrupole relaxation 
tensor is real. Even if time reversal invariance is not 
satisfied, then the he ee, components are real provided 
the site possesses inversion symmetry (Hughes, 1973). 
We shall give the general orientation dependence of Wy 
and Wo for the case of a real relaxation tensor in 
appendix I,since these expressions are somewhat LeEnGEeny. 
Fortunately crystal symmetry will reduce the number of 
independent MBa'B! components (Snyder and Hughes, 1971) 
except for the case of triclinic point group symmetry. 


niieveatciciw lar, for *3 


Na in sodium nitrate where the point 
group symmetry is 3, there are only five independent 


ee atk® components. 


2.7 Mechanisms for quadrupole relaxation 
Se 


We now discuss various mechanisms which have been 
proposed for the quadrupolar interaction of the nuclei 
with the lattice phonons in a solid. The quantities B" 
which occur in eq. (2.26) contain the information 
regarding the lattice vibrations and hence the fluctuat- 


ing field gradients. Expanding BY ina Taylor series 







. 
b= f 2 alced TOR Bie ay ee 
etnsctodmen seabnertakied AL BATES lve een a 
: * f ye Cy use j y \ 7 i] peace ‘h. oes tah ae be 
a oat 
, | 1c. of) sepbattes oo ldon sae ee 


_ se aati a 


4 are Papen 





‘ paeeprat ay ave os od 
a) om ee 


about the equilibrium position, we obtain 


ort 
—s Seren (2039) 
k poke Oe 


where the uy are components of the relative displacements 
and ( 1s Signifies that the quantity is evaluated at the 
equilibrium position, Spin-lattice mechanisms involving 
the term linear in the displacement or quadratic in the 
displacement have been proposed by Van Kranendonk (1954). 
Precxoampleror onesinvolving the linear term is the-so-— 
called direct process in which a lattice phonon is either 
created or annihilated. These phonons have a frequency 
equal to the resonance frequency (2 VON Hz) and, except 
at very low temperatures, constitute only a small frac- 
t10m Of the total phonon population. An example of an 
interaction involving the quadratic term is the indirect 
harmonic Raman process. Here, an incoming phonon 
interacts with the nuclear spin thereby producing a 
phonon of different energy, the energy difference being 
taken up by the nuclear spin. This interaction involves 
phonons from the whole phonon spectrum, and is therefore 
likely to be much more important than the direct process. 
Van Kranendonk and Walker (1967, 1968) have suggested 4 
third possibility which they call the anharmonic Raman 
process. If the lattice potential energy is expanded 


in a Taylor series about the equilibrium configuration, 
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the quadratic term constitutes the harmonic description 
of the lattice whereas higher order terms give rise to 
anharmonic effects. The anharmonic Raman process 
involves the first (cubic) anharmonic term in conjunc-— 
tion with the linear term in eq. (2739), This! inter-— 
action is similar to the harmonic Raman process except 
that an additional virtual phonon is involved. Van 
Kranendonk and Walker (1967, 1968) have shown that for 
both Raman processes, the relaxation rates have the same 


temperature dependence, namely 


WW, = pe] EOme Tea sm0..02 


Vv 
oOo 
e 
ol 


« Te2(7 ~KT*2) for Te (2.40) 


where T* is the so-called reduced temperature T7850 
where oh is the Debye temperature. The constant k is 
very much less than unity. Pietil&d (1968) has shown 
tigi elon Sa in NaCl, the harmonic and anharmonic 
Raman mechanisms are almost the same strength at room 
temperature. Which mechanism, if either, proves to be 
dominant in a particular crystal is therefore still an 


Open question. 
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CHAPTER LiL 


THE = age NMR IN SODIUM NITRATE 


3.1 Crystallographic information on sodium nitrate 
ge as eC OC TUE nL Le Ce: 


We need the structure of the crystal in order to 
predict the electrostatic and magnetic environment of 
the gone Sites.. This explains the NMR of ova in sodium 
nitrate. 

Sodium nitrate belongs to the rhombohedral crystal 
system (space group R3c, Dae and has two units of NaNO, 
per unit cell. As can be seen from Fig. 3.1, the nitro- 
gen atoms are situated at the corners of the unit cell 
and aiso at 1/2, 1/2, 1/2, whereas the sodium atoms are 
ami/4, 174, 1/4 and 3/4, 3/4, 3/4. The nitrate groups 
lie in a plane perpendicular to the three-fold axis, and 
consist of equilateral triangles of oxygen atoms with a 
nitrogen atom at the centroid of the triangle. The 
edges Of the nitrate group triangles are) parallel) to 
the projection of the unit cell edges in a plane per- 
pendicular to the three-fold axis. The orientation of 
the nitrate groups alternate bY fLS 0S Sa long Chis saxis: 

At room Pemparaturer the length (of the body diagonal is 
LO Sze acmewlengin olf nthe unit cells edge is) 53070 Ae, 
and the nitrogen-oxygen bond length is 1.257 A° (Paul 


ancderrvor, £9/1)> “AL room temperature, the motion of 
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The unite Goll of, Sodiumenitra te. mol 
clarity only two groups of oxygen atoms 


are shown. 
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the sodium nuclei is almost completely isotropic, 
whereas the motion of the oxygen atoms is somewhat 
anisotropic (Cherin et al, 1967; Paul and Pryor, 1971). 
Near 276°C, the structure undergoes a second 
erderepiase transition. |This is believed to be an 
order-disorder type of transition in which the orien- 
tation of the nitrate groups becomes disordered. Above 
the phase transition the nitrate groups are statisti- 
cally distributed among the two allowable positions 
OP AleCeiopandyccott, 19/7leand references therein). 
ASwecheresare two sodium atoms per unit.cell in 
the room temperature phase, this means that their 
environments are not identical. Each site possesses 3 
symmetry, the time-averaged field gradient component 
eae being the same at both sites. However, properties 
which involve a fourth-rank tensor are not necessarily 
identical at these sites. The difference between the 
environments of the two sites can be seen in Fig. 3.2. 
The environments are, in each case, similar but are 
differently oriented in space, being symmetrically 


disposed with respect to the three axial glide-planes 


Oo the crystal. 


3.2 The NMR spectrum in monocrystalline sodium nitrate 


Since each soarum nucleus in the unit cell has an 


identical static electric field gradient and since I is 
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Elgunemonc 


The plan view of the nearest neighbour oxygen 
atoms for the A and B sodium sites. The chain 
lines represent the intersections of the glide 
planes with the x-y plane. The DGlnCcipadlwaxcs 
of the M-tensor for the two sites are shown as 


and XprX Also shown is the projection 


MOON Ba 


of the magnetic field Ho onto; the x—-y plane: 
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equal to 3/2, the NMR spectrum consists of a center 
line and one pair of satellite lines. It follows from 
eqs. (2.20) and (2.22) that the resonance frequencies 


are given by 


Vs/¥5 = 1 (0/4) (3 cos*@-1+n sin*6 cos 24) 
De 
ea CaaS Le 260 G3 set) 
£ 2 2 ae 
EAS =e oAG 40s (98 cos 8 — 1) sin 0 (er2) 


where V, is the frequency of the high or low frequency 
satellite, ie is the cénter line frequency and Pp is 
e*Oq/hv.. It is believed that the sodium nuclei are 
situated on a three-fold axis; the asymmetry parameter 
n in eq. (3.1) should therefore be zero. A measure of 
the static quadrupole interaction is e-0q/n, the so- 
called quadrupole coupling constant. At room tempera- 
ture, this has a value of approximately 335 kHz anda 
Cemperaturescockiicientmote-0208 kH2z/eCrm(Andrewletaal, 
1962). At a field of approximately 6700 gauss (the 
strength of our magnet), the unperturbed resonance 
frequency De is 7.55 MHz. Thus, the second order quad- 
rupole shifts are less than 2 kHz and any third order 
shifts (Bersohn, 2952) will be tess than, 10 Hz: 
Inspection,of.eq.4.(3.1) indicates that the second 
opderequedrupole.shait 1si zero ated =.0.0or.90°.. Indéed, 


at all orientations, the separation between the two 


— a 
é ‘to 228ien0 1a ees ete ee tangs - 
i sdztiedves 20 [om eRe. biteit ; 
$ gigit AGS ine 4058) ape 


eA - ~ . oe a] a 
~ 
- 
a w r . 
- . ‘ < } “ ct = 7 
- ~ 7 
- ~- 


tile! co Se all ee 
980, Sad eV — 
| , bape Sat AZ >i _— eee 
say & 7 em 
On ‘ale 
7 
<b avant! igen ts 0) pe ee : 
atuteazesp pi teJa anid) 
: ; 4 LEA, rs UM es}! “bam vetien | 
A) BRS & > ylovsi curt te ater @ fed ins) yam 


fe: do weunadA) 37\ads tt (0-20 serait aeee evaseeapenies 


* | ae lings * . 7 a 
“muy gatiny OOF 3 -didcdamixe=qQae Se) bless fi oA (SPeL 


% 
1. els 


i aurenoze: batgudeegne old , {hedewin smo a cope 





a 
» 


a; i 


Satellite lines is unaffected by the second order shift. 
The relative amplitude of the m< > m-l transi- 


tion is given by (BPP) 


I(I+1l) - m(m-1) . aie, Gill) 


Thus for the case of oe wel obtain 3:74:3 for the relat 
tive amplitude of a satellite, center line, and the 
other satellite respectively. 

We now consider the linewidth of the various 
resonances. The theoretical parameter is the second 
moment Mo, but there is no unique relationship between 
this parameter and the linewidth. HOWGVGr eeLOre ciie 
lineshapes commonly encountered in solids, the following 


relationship is a good approximation (Abragam, 1961): 
Linewidth = 3VM5 (3.4) 


where the linewidth is the peak to peak width of the 
resonance derivative. According to Andrew et al (1962), 
the second moment of the center line for aia in sodium 


nitrate is: 
yi [29.57-2.97 cos 298-9.14 cos 46 + 


+ 6.04(2 sin 20- sin 48)cos 36] 10¢H2- . (45) 


In this case approximately 95% of the second moment is 


from "like" spins; the remaining 5% is caused by the 


an nuclei. Thus, using eq. (2.25), we would expect the 
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second moment of a satellite line to be approximately 


Osso9sofethatwof the .center. line. 


3.3 Orientation dependence of W, and W 


a 2 





The orientation dependence of Wy and Wo fOmsa 
system of identical nuclei at sites possessing 3 
Symmetry is given by (Snyder and Hughes, 1971): 


_ e202 


Ae 
wT BB —{ (My 3394 52. 37Ma333)~ (12M 53 3° 9Mg433) 005 : 


4 
(4M) 111 716M) 3) 348M3333) cos 6 


x * eas) 
16 (M1433 cos 3¢+ Mi123 Sine asin uecos | (Seo. 


TETe od 2 
DE { %yrrat4y 31324559)* (Maya, Giles resie 


4 
+ (My 111741 313723333) COS 8 


# * : ye 
+ 4(Mj 413 cos 36+ Mi 153 sin 3¢)"sin' 6 cos ol. (S37) 


where the z axis coincides with the three-fold axis and 
the. x:and,.y axes, are arbitrary... If the x and y axes 


qa and ee assume 
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L138 


Moe sin 3¢) remains unchanged. The other components 


are rotated, the components M 
different values but the quantity (M cos 346 + 
Mijai° M, 313 and M3333 are unchanged by this transfor- 


mation. It is convenient to choose a coordinate system 
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in which oes is zero; .this means that the x axis 
coincides with a principal axis of the M tensor. 
However, it is impossible to simultaneously satisfy 
this condition for 23Na in NaNO, on account of the 
different environments at the A and B sites. The 
principal axes for these two sites are symmetrically 
disposed relative to the intersection of a glide plane 
with the x-y plane. This is illustrated in Fig. 3.2 
which shows just the nearest neighbour oxygen atoms. 
If we suppose that these atoms vibrate isotropically, 


then from symmetry considerations one would expect 


ands: 


Xar Yq Yp togbesthesxecand vy principalsaxes Hof 


BS 
the M tensor which characterizes each site. Taking 


into account more remote neighbours, x and x 


a’ Ya B, ~B 


are no longer in general the principal axes. However, 
Owing to the short range of the quadrupolar spin-lattice 
teraction (see, for example, Appendix II) the nearest 
neighbour model is likely to be a good approximation: 
for want of evidence to the contrary, we assume that the 
principal axes of the local M tensor do indeed coincide 
Wicu. <4, vn and x 


A 


of the crystal symmetry, we choose the x axis of our 


B’ Yp respectively. To make full use 
coordinate system to coincide with the intersection of 
a glide plane with the x-y plane. Thus, referred to 
this coordinate system, the relaxation probabilities 


for site A will be given by eqs. (3.6) and (3.7) if we 
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make the substitution: 


* x , * 
(M3443 cos 36 + My 453 sin 39) —>Mj 113 cos 3(-$)) 


ie.3) 
where % is the angle between the Xn and x axes as 
SHOW Cel eg. mate. FOrnesiterB, eqs. (3.6) and (357) 


are modified as follows: 


* c * : * 
(M3443 cos 36+Mj 193 Sin 360) > M3443 cos 3(o+60) 2 


(359) 


the A Vand» Bysuciecar Spin systems will be Seong 1 
coupled by spin exchange interactions (Abragam and 
Proctor, 1958; Bloembergen et al, 1959; Andrew and 
Swanson, 1960) due to resonance overlap, since each spin 
System has the same static quadrupole splitting. Also, 
as the heat capacity of each system is identical, the 


observed relaxation probabilities will then be given by 


W, = (W.. + W..)/2 (3410) 
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3 
a0 


a 
2 : ; A 
USs = orm: { 1144! 3394245999) CoM Treetlg gg aucos “o 
+ (M -4M +2M \cos70 + 4M cos 3¢ sin~6 cos 6 
irae 3 Wes 3 353 1113 


Gann) 


where 


Ps * 
a 3 Fey aECOrMs hs Tt 


aBa's! aig Wa 


aremgidentical to those in Wo £orgmuleiphononmrelaxation 


As discussed in section 2.5, the M 


process. Thus the ratio Wo/W, can then be expressed as: 


' ! 1 2 
(144M, 3, 3+2M, 453) +(6-6M3.,) cos Q + | 
+(1-4M! +2M: yeos OF4M. cos 36 sin-6 cos a | 
W, 12153 3333 Tie 
ca ((4+4M! -M! Lh OM. -om' Weos4e | 
V3P3e+3333 13 Be 3g53 55° 
: : 4 4 a 
~ (4-16M, 3, 3+8M,33,) cos 9-16M) 443 cos 3 sin’ 6cos 6 
(3.15) 
where 
i] 
ss Behe a ch ic CA a (3.16) 
' — 
M3333 ~ §3333/@qa11 (3.17) 
' —_— 
NTT Cee (3.18) 


If the quadrupole relaxation is caused. by diffusing 
charges, for example, the Magatg' WILE iot*in general 
be equalin eqs. (3.12yW-and (3.13) (Snyder and Hughes, 


1971). Thus, eq. (3.15) will not be valid in that case. 
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A numerical calculation of the normalized M 


! J 1 
1313’ @3333 and Mj4,3 has been 


given by Reed (1970) and Hughes and Reed (1971) for 


tensor components M 


a nearest neighbour point-charge model of sodium nitrate. 
The detailed expressions are somewhat lengthy and are 


therefore given in Appendix II. 


3.4 Spin-lattice relaxation transitions 


The oone quadrupole coupling constant is suffi- 
ciently large that the three resonances are well resolved 
at most crystal orientations. For an accurate calcula- 


Gon Of the relaxation behaviour of the 23 


Na spin system, 
it is necessary to take into account the quadrupole 
shifts of the energy levels. However, for large magnetic 
fields, the ecensaorde: quadrupole shifts are small 


compared with the first order shifts. Thus, we can write 


the resonance frequencies as 


Ae vo (i+ 9) (3.20) 
aga s (3 i218) 
vata, (lysed) Bt (3-22) 
Since I = 3/2 for ater the quadrupole relaxation 


probabilities are given by (see eqs. (2.31) and (2.32)) 
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Wa, = 0 (3.24) 
Woeme  a a oS (m2 sy 


However, for a population dynamics Calor vor Us s 
necessary to bear in mind that the upward and downward 
transition probabilities are not exactly equal, since 
they are related by the Boltzmann factor. We therefore 
choose to call the upward transition probabilities Wy 


and Wo, and the corresponding downward probabilities 


W, (1 + A+é6é) and W, (1+ 20+ 6) where 
a= hv o/kT, A C3re2 6) 


These transitions are shown in Ub eik) Sikes 

In general, one may also expect some magnetic 
relaxation, due, for example, to paramagnetic impurities. 
The upward transition probability for this case is given 


by (Andrew and Swanson, 1960): 


Meer = W., (I+m) (I-m+1) ia eeehy) 


and the downward transition probability is obtained by 
multiplying eq. (3.27) by the appropriate Boltzmann 
factor. These transitions are also shown in Fig. 3.3. 
If the resonances are not well resolved, 
additional transitions will occur due to the so-called 
Spin-exchange interaction (Abragam and PVOCtOr, «1958" 


Bloembergen et al, 1959; Andrew and Swanson, 1960). 
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The spin-lattice relaxation transitions 
for the case I = 3/2, where the quadrupole 
Splvtceing is subbiclen pacOmre sO! Dem crete 


exchange transitions. 
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Since these transitions introduce an additional coupl- 
ing between the energy levels, they make it more 


difficult to interpret spin-lattice relaxation measure- 


ments. 
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CHAPTER IV 


EXPERIMENTAL METHODS FOR STUDYING THE 


QUADRUPOLE SPIN-LATTICE RELAXATION 


4.1 Steady-state double resonance method for measuring 
WO 

It was first demonstrated by Pound (1950) that 
it is possible to distinguish between Magnetic and 
quadrupole relaxation mechanisms by means of a double 
resonance experiment on a quadrupole-split spectrum. 
For the case I equals 3/2, Pound showed that if Wy equals 
W5,and W3 equals zero,the population difference of a 
Satellite (i.e. its amplitude) is increased by 50% when 
the center line is severely saturated. For the case 
where magnetic relaxation is dominant, Pound showed that 
Saturating one resonance would cause no change in the 
others. Later, Abragam (1961) generalized the theory 


and showed that the satellite enhancement E is given by 
BE = (1+ 2x)/(1+ x) (4.1) 


where x is Wo/Wy- Thus, a measurement of the satellite 
enhancement should provide a value of Wo/W) - tt Ave 
obviously essential that nuclei involved must be at 
sites of lower than’ cubic symmetry.9 It°is also, highly 


desirable that the sample should be a single crystal, 
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SO as to prevent smearing of the satellites. Abragam 
showed that Wo/Wy could also be obtained by severely 
Saturating a satellite and measuring the amplitude of 
the center line or the other satellite, the enhance- 


ments being respectively equal to 


ea] 
I 


(3 + 2x) /(2+x) (42) 


and 


3] 
I 


OP 2at eX), (4.3) 


We now present a simple physical picture which 
the author has found useful in explaining the enhance- 
ment of the satellite when the center line is severely 
Saturated. The spin-lattice relaxation transitions 
tend to keep the spin system at the lattice temperature, 
whereas an rf field applied on resonance tends to 
decrease the population difference of the two energy 
levels involved. When a large saturation field is 
applied to the center line, its population difference 
rapidly decays to a smali value (zero in the case of 
complete saturation). It can be seen from Fig. 3.3 
that the population difference of each satellite line 
must correspondingly increase rapidly to 1.5 Noe where 
Ny is the equilibrium population difference in the 
absence of applied rf fields (in this discussion, we 
neglect the small correction factor é).... The Wy 


mechanism will tend to reduce the satellite population 
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difference to Noe whereas the W. mechanism will tend 

to increase it to 2n.- This last point can be seen by 
remembering that the W, mechanism couples levels 
separated by Am= +2 (see Fig. 3.3) and that the center 
line population difference is maintained at zero by 

the saturation field. Thus, the pee Sten Tas ate 
population difference of the satellite line depends 
upon the relative strength of Wo and Wi: When W/W, is 
very much greater than unity, the population difference 
approaches 2nor anol. L W/W, is very much less than unity 
the population difference approaches Noi these values 
can be compared with the predictions of eq. (4.1). 

To obtain accurate W/W, values from enhancement 
measurements, the observing power of the spectrometer 
must not significantly perturb the spin system. The 
effect of a finite observing power is to make Wy "less 
effective". Thus, the measured enhancement will be 
larger than the true value. Theoretical expressions 
for the enhancement in terms of W,/Wy were obtained by 
Hughes (1966) who also showed that the perturbing effect 
of a small observing power could be allowed for by a 
simple linear extrapolation to zero observing power. 
Reed (1970) showed experimentally that the enhancement 
did, indeed, vary linearly with observing power over 

23 


quite an extended range for the Na spin system in 


sodium nitrate. 
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The method for calculating the enhancement 
(Pound, 1950; Abragam, 1961; Hughes, 1966; Reed, 1970) 
consists in applying the principle of detailed balanc- 
ing to the spin system, in order to obtain rate equations 
for the energy level populations. We shall not give 
the details of the steps but instead quote the steady 
state enhancement expressions given by Reed (1970). 

The first case (Case I) is where the high fre- 
quency satellite is observed with spectrometer power 
Eee and the center line is completely saturated: the 


enhancement is given by 


_ — (1+2x) 7 6X ‘ 3xXY ee fobs (4.4) 
i Tees (1+2x) eo a) 


where y is W3/W, , the normalized magnetic spin-lattice 
transition probability. 

Case II is where the low frequency satellite is 
observed and the center line is completely saturated. 


In this case, the enhancement is 


Shag? 
Siti 2x) ax bye 3xy obs 
Stn Taney fa (tax) - (iexhGlezx)) | w2lac ae e 


Case III is where the center line is observed 
and the high frequency satellite is completely saturated. 
itGelsefoundsthat: 


E he (3+2x) ioe 6 (1+x) : aren rt 2 (l+x)P og 
EL tee (2-3) (3+2x) x(24x) (342) “ees cae a 


(4.6) 
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Case IV is where the center line is observed and 
the low frequency satellite is completely saturated. 


Here one finds that: 


peer eee y ape) ope Noe 
Te (2 +x) (34+2x) x(2+x) (3+2x) (2+x) ei r 


The expressions for re co Esty avee Val Veron any 
Mea BUC Olax Out yen Oo. and Dae are small correction factors 
and, thus, must be much less than unity. These expres- 
Sions also assume that the appropriate resonance is 
completely saturated. This obviously cannot be achieved 
with a finite saturation field. Thus, we should correct 
the enhancement values to take into account the incom- 
plete’ saturation: .When 6, Rage and syeanehequals tol zero, 
the enhancement as a function of saturation power Pra 
is given by: 
2 422P eee x) yc 
ae ees 2) ad ae 
1+2P_.,[(1tx) /x] 
POLBCASCS Ee and sil. ihesexpression for Cases) llimanawi. 
is similarly given by: 
1+ 2P_6,[(3t2x) /x] 


E ee 
PRET AN Be 
1+ 2Pli,[(2tx) /x] 


(4.9) 


The validity of the equations (4.8) and (4.9) was con- 
firmed by Reed (1970).. Thus, they can be used to make 
small corrections to the enhancement for the fact that 


the saturation power is not infinite. 
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Since the expressions (4.4) to (4.7) contain 
two unknown quantities, x and y (6 can easily be 
obtained from the satellite splitting) hit is necessary 
to make two independent enhancement measurements. 
From the form of the equations, it can. be seen that 
one must be either Case I or II, whereas the other 
RUStspet case Lior et Vow ele Veeindeed!) fortunatesthat 
Wecieers drrrerentiy inseqs ©44.4)Wand (495)? than it 
Re in eqs. (4.6) and (4.7); it would otherwise be 


impossible to determine x and y separately. 


4,2 Transient double resonance method for measuring Wy 





The steady state enhancement experiment just 
described gives only the ratio W./W, - However, it is 
also useful to know the magnitude of Wy (eke We separately. 
Puémencoret ical expressions (eqga (3.12) and (3.13)) 
give the orientation dependence of these quantities, 


and it is therefore only necessary to determine W, or 


aa 
Wo at one crystal orientation. 

Ingprinciple it is@possible  tovobtain Wy from a 
steady state progressive saturation experiment (Hughes, 
L665 However, the accuracy of such an experiment is 
very limited and a direct measurement of the time 
dependence of the spin system is more fruitful.- This 


may be achieved by producing a non-equilibrium state of 


the spin system and observing its recovery. 
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Expressions for the transient behaviour of the 
Spin system have been given by Niemel& (1968) for the 
approximation S orP re = 0. The transients are in 
general multi~exponential functions of time, and the 
task of determining the individual decay constants 
from experimentally obtained data is One Of considerable 
complexity (Smith and Cohn-Sfetcu, 1973). However, 
using the expressions given by Niemeld (96S) ert acan 
be shown that the recovery of a satellite line after 


removal of the saturation power from the center line 


should be of the forn: 


A+B exp =At |. (4.10) 


We carried out a detailed Calculation of this case, tak- 
ing into account first order 6 and ae corrections but 
assuming that y is negligible. The population difference 
n, corresponding to the high frequency satellite has the 


following time dependence: 





3 (2+x)P 
a obs 
n, =n (1+) h 2 (1+X) | 
3 (x2+x-1)P 
x ODS 
ae 2 (14x) 2 Jexee~(amy+3P9, gm) 


obs 
y 


+ 





fees 
n exp{-(2W,+2W,+3P_, w.)t} . (4.11) 
a) 2 (14%) iL 2 obs 1 


We note that there are now two exponential terms rather 
than one. However, the amplitude of the second expon- 


ential is proportional to ade and is therefore small, 
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provided one uses a low observing power. We also note 
that the decay constant of the dominant exponential is 
a linear function of Pobs* Thus, by measuring the decay 


constant as a function of Po Onensnouldebe vableasto 


bse 


determine Wy by use of a linear extrapolation to Bebe 
equals zero. The effect of magnetic relaxation was 
maker ynco account for the case ee = 0. It was found 
that the decay constant of the dominant exponential is 


equalaco: 
W, (l+ x+7y+ | (1-x+5y)* + 4y (x-1)) (4a) 


rather than 2Wy as before. 

Thus, we see that Ww, can be measured by monitoring 
the recovery of the satellite amplitude after removing 
the saturating power from the center line. The decay 
constant is then plotted asta function of "spectrometer 
power and the intercept is equated to expression (4.12). 


TE might be thought that W, could also be obtained 


al 
from the transient behaviour of a satellite when the 
center line is first saturated. We shall now describe 


why this is not practical. -The satellite population 


difference has the following time dependence: 


Nees 
Tes] +n, [gt] expt- (2how,t) 


Oo 
= “= exp{-6P__,W,t} ° (4.513) 


aul 

















ra. iy fy vio Sly ney iSIhtid BPE , tape *fasS 
see : ¢ 
ic Si omideh sae tee Ste. eee, a fais | qi 
wes on el eae eee 
BASC 4 * 


pe ife off )20h See eo Rea eae Seaenad| 


17. 39° Setende qete saa aa 


ted feut 


—— ee I [ .* a * C): te - 


“! cA pis nats 2a 
rh 66 es 1 Jac ape oes, ope? =~ 
ws silt; 2c ~wreavaoet is", 
ey tall Avion 4 i$ eae gots sraiviea ote 
ae we Si ley eee ai 209:12009) 
t 73 Oo! (beynvpy Se Doane ate —_— : vem : 
bor <3 4 Wiad > shns sree aa: val dea #2 
att obit Leta: p. to. 2eelaetog greens Pak 


Re sees es 


= 
Gg 


i -s he 


ie 


a _ ' 
7 = 
ay 


D2 


The last term in eq. (4.13) represents the rapid rise 
COR. No described in section 4.1. Observation of the 
time dependent behaviour corresponding to the second 
term would give the value of (2+x)W, . However, the 
amplitude of this term will be small unless x is very 


different from unity. 
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CHAPTER V 
APPARATUS AND EXPERIMENTAL DETAILS 


5.1 Measurement of the NMR enhancement 
a ee es CH TICeMel eC 


The resonance enhancement E is determined by 
measuring, as a function of spectrometer observing 
power Bee the increase in amplitude of one component 
of the three-line spectrum when another component is 
severely saturated. The perturbing effect of the 
observing power is accounted for ee a linear extrapola- 


tion procedure of the E-Po curve to Ps = 0 (Hughes, 


bs bs 


1966; Reed, 1970). 

For a continuous wave experiment, the experimen—_ 
tal constraints on the spectrometer are, firstly, it 
should be variable frequency-fixed field rather than 
vice versa, since a variable field implies that the 
saturation frequency would have to track the resonance 
frequency. Secondly, it should be a single coil type, 
because of the need for an orthogonal coil to provide 
the saturation field. Thirdly, the spectrometer should 
Operate at low power levels. such that the transition 
rate in the spin system caused by the spectrometer is 
small compared with the spin-lattice relaxation rate. 

The first constraint. implies that Rollin (1946) 
and bridge spectrometers (BPP) are unsuitable, while the 


second constraint rules out the crossed-coil variety 


Sp) 
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(Bloch, Hansen and Packard, 1946). The marginal 
(Watkins and Pound, 1951) and Robinson (1959) self- 
oscillating spectrometers satisfy the three constraints. 
The Robinson type was chosen on account of its stable 
and reliable operation at the low power levels required. 
Moreover, it responds essentially to only the absorp- 
tCtonesignaly({Hughessand.fmith,.19/1)-.. 

Reed (1970) measured the peak-to-peak amplitude 
of the resonance derivative by conventional lock-in 
detection. and sinusoidal field modulation in conjunc- 
tion with a linear frequency sweep through the resonance. 
Reed's method is inefficient, since only a small frac- 
tion of the time is spent near the resonance extrema. 
While this disadvantage may be overcome by measuring the 
resonance amplitude only at an extremum, there remains 
the problem of maintaining the appropriate spectrometer 
frequency. 

An alternative method is to use two lock-in 
detectors, one operating at the modulation frequency 


Wea to lock the spectrometer onto the resonance, the 


54 


other operating at 20 to measure the resonance amplitude. 


To maintain the spectrometer frequency, the error signal 


should be zero at the desired locking point and should 


a 


Lock-in detectors are also known as phase-sensitive 


detectors, correlators or synchronous demodulators. 
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also be an odd (preferably linear) function of frequency 
deviation. These conditions are fulfilled by the out- 
put of the On lock-in detector at the center of symmetri- 
cal resonances, for all modulation amplitudes. For good 
signal-to-noise ratio, the spectrometer frequency should 
be chosen to give a maximum output from the 20 lock-in 
detector. This condition is again fulfilled at the 
resonance center. In fact, the resonance center is the 
only position having all the desired properties and 
still allowing freedom to optimize the modulation ampli- 
tude. 

In the limit of small modulation amplitudes the 
On and 20, lock-in detector signals are proportional to 
the first and second derivatives respectively of the 
lineshape (Wilson, 1963); these are shown for a Gaussian 
lineshape in Fig. 5.1. 

Most of the resonance amplitude values reported 
in this thesis were obtained using this arrangement in 
which the 20 lock-in detector output was measured while 
the Wn lock-in detector maintained the spectrometer 
frequency at the resonance center. The enhancement 
values are equal to the ratio of the resonance ampli- 
tudes in the presence and absence of the saturation 
power. The error due to any de offset in the circuitry 
following the lock-in detector was eliminated by perio- 


dically reversing the phase of the 20, reference signal. 
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Fiquremoms 


A - Gaussian lineshape. 

B - First derivative of the Gaussian 
lineshape. 

C - Second derivative of the Gaussian 


lineshape. 
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The measurement procedure is illustrated in Fig. 5.2. 
This ten-minute"saturation-cycle"was repeated until 
the desired precision in E was obtained. The enhance- 
ment was measured at a number of fiprecent values of 


Pp so that the value of E at Eee equals zero could 


obs 
be determined by the linear extrapolation procedure. 

The small correction factor to take into account 
the non-infinite saturation power was determined by the 
following procedure. We measured the enhancement for 
three or four values of saturation power (at a low 
value of aa ey and the value of E at infinite power 
was determined using equation (4.8). This value was 
typically 0.15% greater than the value obtained at the 
finite saturation power used in the previous measure- 
ments. We then used this to correct the extrapolated 
value of E at oe = 0. 

For most of this work, the Robinson spectrometer 
was locked to the high frequency satellite and the 
Saturation power was applied to the center line. 
Approximately forty-eight to sixty hours of data were 
collected at each crystal orientation in order to attain 
the desired precision in E. Thus, the advantage of a 
mini-computer controlled system operating twenty-four 
hours per day is obvious. 


To obtain accurate and reliable values of the 


enhancement, we require the spectrometer frequency to 
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OUuEDpUL OLSthe 20 lock-in detector for 
one complete "Saturation cycle". The 


bandwidth of the low-pass filter is 1 Hz. 
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be locked, the magnetic field to remain constant and 
the spectrometer sensitivity to NMR to be unchanged 
throughout the saturation cycle. A large amount of 
ancillary equipment was required to achieve these 
goals, as will be described in the sections that 


follow. 


9.2 Description of the apparatus 
Seow Block diadram 


The block-diagram of the system is shown in 
Fig. 5.3. Excluding commercial equipment, the system 
contains approximately 600 integrated circuits, both 
digital and analogue, plus numerous discrete components, 
semiconductor devices and thermionic valves. The indi- 
vidual circuit diagrams of the blocks will therefore 
not be given. Instead, we shall describe their function, 
operation and important features. 

Most of the digital and analogue circuitry was 
assembled on 4%" x 6%" circuit cards which are plugged 
into multiple -pin edge connectors, thereby assuring 
ease of access for maintenance. The construction and 
development of the system extended over a period in 
excess of three years and included a considerable amount 
of wedesign to obtain tthe final working version. 

The central feature of the system is the (System) 


Buss. LES LUnCE LOM LS to enable the various devices to 
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Block diagram of the minicomputer-controlled 
double resonance spectrometer system. The 
arrows indicate the direction of information 


or signal flow. 
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communicate data and commands to each other. It con- 
SiptosOlelondatamiines we loamonitonue lines, aS well as 
several address and strobe lines. Each device is 
eonneccted eithersdirectly or indirectly to the Buss. 

The interface is a bi-directional device enabling 
information to be transmitted between the Buss and the 
computer. As an example of the operation of the system, 
we shall consider one method of frequency counter con- 
trol. The computer sends a command on the address and 
strobe lines to the frequency counter to start counting. 
One second later, when the frequency counter has com- 
pleted its measurement, an interrupt signal is sent to 
ENewcompuUcer, whrchamay Or May Not act upon it. “If it 
does, it may command the frequency counter to transfer 
the digital representation of the measured frequency 
onto the monitor lines and thus into the computer. 

The blocks are described in detail in sections 
5.2.2 to 5.2.17. Meanwhile, we describe briefly their 
functional relationships. 

The NaNO, Sample is positioned within the spec- 
trometer coil, which in turn is located in the aluminium 
probe. The mp NMR signal is detected by the spectro- 
meter in conjunction with the On and 204 lock-in 
detectors. The "modulation synthesizer and phase control" 


provides the signal feeding the modulation coils, and 
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the reference signals for the lock-in detectors. The 
phase of the reference signals relative to the modu- 

lation signal can be controlled by the system or can 

be adjusted manually. 

Bache lock-inidetector, output,.is, fed, to. the 
analogue multiplexer,and then digitized by the analogue 
to digital (A/D) converter before being coupled to the 
Buss. 

The spectrometer frequency is controlled in the 
following way: the frequency converter and counter 
together provide a digital representation of the 
spectrometer frequency. This information is compared 
with the "correct frequency" stored in the computer, 
and any correction is fed to the frequency control D/A 
converter. The analogue output biases a variable 
Capacitance diode in the oscillator tuned circuit. 

The rf power source drives a pair of coils 
within the probe which provide the saturation field. 
Although these coils are nominally orthogonal to the 
spectrometer coil, some of the saturation signal 
unavoidably leaks into the spectrometer. We have 
shiown| theoretically (Hughes etal ,elo7 4). thatsthe 
sensitivity to NMR of a Robinson spectrometer should 
be unaffected by a small leakage. The sensitivity of 
our spectrometer, however, was found experimentally to 


vary quadratically with the injection voltage; this 
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is attributed to non-ideality of the spectrometer. In 
order to overcome this problem, the leakage is largely 
eliminated by the flux balance network, which is based 
on the design of Hughes and Reed (1970). The deviation 
from perfect balance is monitored by the misbalance 
detector which, on command, can inform the computer of 
Teotmagnitude.- this information is then used by.the 
computer to readjust the balance network. 

The NMR field-lock spectrometer maintains the 
magnetic field in the magnet gap at a predetermined 
value throughout the course of the experiment. 

The manual operation console enables the operator 
to interact directly with the system, thus by-passing 
PheecOMbuULer mwalCheis the real-time control center of 
the system. The teletype provides a means of entering 
programs into the computer as well as printing the 


experimental results. 


Seek Probe 


This consists of an aluminium body with nylon 
inserts to provide mechanical support for the spectro- 


meter and saturation coils. These were wound of silver 
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wire rather, than-copper;, since, the Cu NMR signal 


occurs close to the 23ua resonance. To increase the 


"NMR filling factor" the spectrometer coil was glued 
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to the inside rather than the outside of a nylon holder. 
The position of the connecting leads was adjusted so 

as to minimize the net coupling with the saturation 
Coils. The size.and inductance of the latter were 
chosen to maximize the saturation field produced at 

the sample. The spectrometer, probe and balance net- 
work were interconnected by triaxial cable, as shown 

in Fig. 5.4, so that the rf return currents would not 
flow in the outer shield. 

The modulation coils consisted of a hundred turns 
of copper wire, which were wound on nylon formers 
outside the aluminium probe. The leads to the modula- 
tion coils were also kept outside the probe body to 
minimize coupling with the sensitive rf CEECUmLT yy. 

The sample was mounted in a nylon holder which 
enabled the crystal to be rotated about two mutually 
perpendicular axes, as will be described in section 


sie Ue 


5.2.3 Robinson spectrometer 


A Robinson type NMR spectrometer (Robinson, 1959) 
consists of a high-Q parallel-tuned L-C circuit followed 
by a low-noise rf amplifier, limiter and feedback net- 
work to sustain oscillations. Changes in oscillation 


level caused by NMR absorption are monitored by a 
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Figure 5.4 


Schematic diagram of the wiring 
associated with the NMR probe and 


the Balance Network. 
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diode detector followed by a low noise audio frequency 
amplifier. 

Our design is a modification of that given by 
Howling (1966) and the essential features are shown 
Pel ig. S.o. limiting 1S achteved in two stages using 
Schottky diodes. These were used because of their low 
noise, low capacitance, and sharp "turn-on" charac- 
teristics. Five switchable resistive dividers enabled 
=the oscillation level to be set over the range 5 to 
50 mV RMS. The feedback network comprises ten 20 Kf 
resistors connected in series and located in the center 
of a large metal box. This arrangement provides a iow- 
loss feedback path which is essentially resistive. 

Ground loops were avoided by having only one dc 
connection between the system and the spectrometer, 
this being the ground connection to the variable capa- 
Citance diode. Other dc connections were avoided by 
using transformer coupling for the audio output and 
optical couplers for the misbalance detector and fre- 
quency converter outputs. 

A major problem was the spurious mixing of rf 
and audio frequencies. We note Gaels if the oscilla- 
tion signal is modulated at Ww, OF 2u ar these will 
appear as fake outputs at the Wn and 20 lock-in 


getectors., The latter is particularly troublesome, 
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Figure 5.5 


Block diagram showing the essential features 
of the Robinson spectrometer used to measure 


NMR amplitudes. The sample is situated within 


theasinaductorm in @). 
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Since it introduces a systematic error in the deter- 
mination of enhancement values, as discussed in section 
5.4.3. These effects were ameliorated by extensive rf 
£vUlteringtand burfiering ofCall#connections tothe 
Spectrometer. 

Some time was spent modifying the spectrometer 
so that adequate NMR signal-to-noise ratios were ob- 
tained! —In “particular, the gain ‘of rthe*rf amplifiers 
was made larger than one hundred so that noise origin- 
ating in the detector and audio stages is unimportant. 
The signal-to-noise ratio is proportional to the square 
root of the in-situ Q factor of the sample tuned circuit. 
We found that this” O-factor can be significantly 
reduced by seemingly innocuous components such as the 
502 grid-stopper used by Howling (1966) in the input 
stage of his spectrometer. The effective NMR filling 
factor was improved by minimizing the inductance of the 
connection leads to the sample coil. The triaxial shield- 
ing arrangement mentioned in section 5.2.2 also served 
to reduce the sensitivity of the spectrometer to exter- 


nal interference. 


5.2.4 Spectrometer level monitor 


A "three-and-a-half-digit" dc voltmeter measures 
the spectrometer detector voltage, which in turn is 


related to the spectrometer level. This voltmeter was 
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calibrated in the following way: The Marconi oscillator 
(see section 5.2.5) was connected across the spectro- 
meter tuned circuit and the voltmeter was read as a 
Function» of the Marconi oscillator output. “Thus; the 
spectrometer level is the equivalent Marconi output which 
would give the observed detector level. Since the 
characteristics of the spectrometer may be frequency 
dependent, this calibration procedure was performed 
every time a new crystal orientation was selected. 

The output attenuator of the Marconi oscillator 
was in turn calibrated by a "true RMS thermal converter" 


which measures rf voltages to better than 1%. 


5.2.5 The rf power source 


Preliminary measurements using a Takeda-Riken 
Model 2150 frequency synthesizer showed that it was 
too noisy to use as the source of the saturation signal; 
this was because the balance condition of the balance 
network was frequency Gesendent. Thus, even though 
there may be no coupling between the saturation and 
spectrometer coils at the saturation frequency, this 
is not the case at other frequencies. In particular, 
noise generated by the saturation signal source near 
the spectrometer frequency will be coupled to the 
spectrometer. A Marconi TF2002 L-C signal generator 


was found to be sufficiently noise-free. However, it 
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lacked the required frequency stability. A Marconi 
TF2170B digital synchronizer was therefore used to 

lock the signal generator, thereby achieving a frequency 
Stability of +1 Hz. The Marconi output is then fed, 

via the system-controlled co-axial on/off switch to 

an E.N.I. Model 310L ten watt rf power amplifier. To 
provide electrical isolation,the logic signal to the 
switch is passed through an optical coupler. The power 
amplifier output is then fed to the saturation coils and 


the balance network. 


5.2.6 Balance Network 


Although the connection leads of the spectrometer 
coil within the probe were adjusted to minimize the 
coupling with the saturation coils, a non-zero leakage 
remained. This was found to be approximately orthogonal 
to the saturation voltage; we attribute this to the 
eddy currents induced in the aluminium probe case. 

The flux balance network is based on a design 
of Hughes and Reed (1970), and operates by synthesizing 
a cancelling leakage from in-phase and quadrature mode 
signals, each of adjustable amplitude. All the leakage 
signals are coherent, since they are derived from the 
same rf power source. As shown in Fig. 5.4, injection 


signals are coupled to loops in the spectrometer circuit. 
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Reversible synchronous motors and a series of reduction 
gears rotate grounded aluminium vanes between the loops, 
thereby adjusting the amplitude of the couplings. The 
two modes were set precisely orthogonal to one another 
by adjustment of a variable delay line in the orthogonal 
mode circuitry. This arrangement gave two orthogonal 
modes with little interaction between them, thereby 
enabling the balance condition to be reached in a small 
number of steps. 

Because of the imperfections in the spectrometer, 
mentioned in section 5.2.3, a net leakage of less than 
Peart lin 107 of the saturation voltage was required. 
This placed great demands on the stability of the 
balance network. Consequently, all loops, vanes and 
Wir tiCmewereuTi gi 0ly constructed CO minimize micro-— 
phonics. One major source of instability was Joule 
heating of the delay line. This was largely overcome 
by spn arbi the coupling, thereby reducing the current 
flow through the delay line. 

Amiepaitat we” correction to a mode corresponds 
to a vane rotation of only three minutes of arc. High 
quality bearings and anti~backlash gears were therefore 
used to ensure smooth and reliable rotation. Even so, 
backlash still remained a problem, but was solved at 


the software rather than the hardware level. 


Tes) 





7 being oa ve 






















7 ae 

¢ bas Bhs gone LD ye sual ‘7 
jie En Bebe ore otnton ame 
sty J bike’ a3 “patsy fhe vdotadd 
stay aohowiond & 


| .. 
st & Lo retisenthe Bay 


so yards Vleg tees See- 


aie el 


r 
. Pi - 
aoe ifs Le SF) 


4s igouske 
a “< 7 


atsoliods iw ¢sbam 


4 od) hety Lomo anue ed waft paris 
> ; a eget In: 
teqeb ete Be javesa 
Ts :.9,@ doldooe at Bengise 
: ee to ‘on ni Haq £ 
¢ ehopbe> tame Saget 
® ti tooo peRiee  azowtem Some 
400 eibiols e168. 
tia .4 ars Ot wut 280 
cow - OLR wett yeleb odd 16 pi 





stenmy ,vakiques ait? ‘pale. 


ones sgl ed. a 
heyy a cha 








76 


The sense and the duration of the rotation of 
each motor was controlled by a bank of relays, which 
in turn were activated by the system. Optical couplers 
between the system and the relays provided electrical 


isolation,and minimized interference effects. 


5.2.7: Misbalance detector 


Any residual saturation signal entering the spec- 
trometer will be mixed with the spectrometer signals by 
the nonlinear elements, thus producing, sum and difference 
frequencies. The latter component, which is less than 
170 kHz in our experiment, is amplified, filtered, 
rectified and smoothed in the misbalance detector unit. 
To avoid the dc loop created by connecting the output 
of the misbalance detector directly to the analogue 
multiplexer, we used an ac method of coupling into the 
system. The slowly varying dc voltage modulates the 
pulse width of a 100 Hz square wave carrier, which is 
then passed to an optical coupler and hence to the 
system. This pulse train is used to gate a 10 kHz 
clock pulse. By counting the number of clock pulses 
‘per Carrier cycle, we obtain aybinatry number Tcrepresen= 
ting the orioinaledcw voltages touaiwaccuracy, On (le. 

An interrupt signal is sent to the computer whenever 
this number is greater than a predetermined reference 


number. 
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Previous work in this laboratory has shown 
(Hughes et al, 1973) that the precise relationship 
between the injection voltage and the change in 
spectrometer NMR sensitivity is a complicated function 
of spectrometer circuit parameters and the difference 
frequency. The misbalance detector was, therefore, 
calibrated by the following procedure. The balance 
network was misbalanced by trial and error until the 
NMR spectrometer sensitivity changed by 1%. Then, 
assuming a quadratic dependence of NMR sensitivity 
change on the injection voltage as found by Hughes et 
pists) mcs mjecelong voltage was, reduced to the 
level which would give a sensitivity change of 0.05%. 
The 1% sensitivity change was chosen to ensure the 
wolidutysotthe extrapolation. The output, of the mis— 
balance detector was noted and the experiment repeated 
as a function of difference frequency and spectrometer 
levers tiusm enablingathe interrupt, level to be set to 
a level such that the change in the NMR sensitivity is 
no more than 0.05%. The good signal-to-noise ratio of 
the Pees resonance from a sample of metallic sodium 
facilitated. the calibration procedure. 4 The computer 
program, written for the enhancement experiment on 


NaNO was used without modification for the measure- 


Bilt 


ment of the change of NMR sensitivity. 
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5.2.8 Frequency converter 


The frequency counter described in the next 
section consists of a Series of divide-by-two stages. 
The outputs of these stages will therefore, in general, 
have a component at the spectrometer frequency. More- 
Over, the amplitude of this component depends on whether 
or not the frequency counter is actuabiy counting. @Any 
radiation from the frequency counter entering the 
spectrometer will therefore appear as an additional 
and variable feedback path. 

To avoid this problem, the spectrometer frequency 
was heterodyned with a crystal-controlled OSseil lator, 
whose frequency was 9.76 MHz. This local oscillator 
frequency was selected to ensure that no harmonic of 
the intermediate frequency (approximately 2.2 MHz) could 
ever coincide with the spectrometer frequency over its 
Operating range. The output of the frequency "converter 
was fed to the frequency counter via an optical coupler 
providing electrical isolation. 

To reduce radiation at the spectrometer fre- 
quency, the frequency converter was mounted in a steel 
box near the spectrometer, so reducing lead lengths; 
this unit also had multistage rf filtering, or the, dc 


Supply lines and in the ac power line to it. 
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5.2.9 Frequency counter 


The system-controlled frequency counter has a 
resolution of 1 Hz and a count time of one second 
derived from a stable 1 MHz crystal oscillator (/RMS 
Model M76BR). After the counter has completed its 
measurement, it may be restarted within 1 millisecond, 
unlike most commercial counters in which a wait of one 
second would be required. The count number is in binary 
format since this is a convenient representation for 
the computer to manipulate. (The binary coded decimal 
B.C.D. format is convenient for a visual display but 
not so convenient for computer manipulation.) The 22 
bit binary equivalent of the intermediate frequency is 
too long to be handled by the computer software, there- 
foresonlysthe-least) significant 12 bits. are used.-iThe 
resulting ambiguity in frequency is not a problem 
provided the spectrometer frequency changes are less 
Gnane!t bits/sec (20480iz/sec) - In practice this con— 
dition is easily satisfied since the spectrometer 


instability is seldom in excess of 10 Hz/sec. 


eeel0) Frequency sadjustmentsclrcuicry 


In a semiconductor diode, the width of the deple- 
tion layer, and hence the capacitance, is a function of 


the bias; thissproperty. is,utilized in variable 
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Capacitance diodes. If such a device forms part of 
thie spectromecer input tuned circuit, it provides an 
electrical means of frequency control. For this par- 
BYCULaAlL@dapDLECatLol,. the OUtTpUL Of a 12 bit D/A 
converter reverse biases the diode. The bit size 
corresponds to a frequency change of 2 HZ, giving a 
total range of just over 8 kHz. Since noise voltages 
applied to the variable capacitance diode frequency 
modulate the spectrometer, the control voltage was 
passed through a multistage 10 Hz low-pass filter, 
thereby reducing any spurious frequency modulation, at 
ii and 20 ins particular. 

Data from the computer is fed into the frequency 
control device via the data lines, and then stored in 
12 data latches. Their outputs are connected to the 
D/A converter and their contents may be ascertained by 


the monitor lines. 


See Leet ie Wn and Zu, lock-in detectors 


Lock-in detectors are commonly used in various 
branches of spectroscopy where the phenomenon under 
investigation can be modulated by a periodic signal. 

In general, the spectrometer output can then be written 


as 
Mie = L a cos (nut + >.) (5b) 





where w is the modulation frequency. If Nie is multi- 


plied by a reference signal 


Vo D B, cos(nwt +p.) (S020 


of the same periodicity, the dc component of the result- 


ing signal is 


NS) 


» A,B, Cos (¢,- v,,) : (Sao) 


This process is essentially the determination of the 
cross-correlation coefficient of We and Mt tae frou 
components of Ve except the m'th are rejected by 
Suitable filtering before the multiplication stage, 


the dc component is 


l 
5 avon cos(¢. vn é (524) 
This has a maximum when Coreg —sUmOrmrlCUl — MOLeCOVeL, 


hg & Ba is constant, this output is a measure of the 
amplitude of the m'th harmonic of Vo: We note that 
for the latter case, the shape of the reference signal 
is unimportant, as long as the amplitude of the m'th 
harmonic remains constant. 

Our lock-in detectors operate on this principle, 
and the design has been described in some detail by 
Boyd) (19/5). Each detector Consists vol an inpuc stage, 


a band-pass filter and amplifier, a correlator anda 
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post-correlator low-pass filter. 

The choice of modulation frequency was governed 
by the following constraints. It should be high enough 
POsovolds the flicker noise” or 1/f noise. region. 
However, it should be small compared with the NMR 
linewidth expressed in frequency units (Haworth and 
Richards, 1972). Also, frequencies near 60 Hz and 120 Hz 
should be avoided. Indeed, it is desirable that har- 
monics of the modulation frequency should not coincide 
with harmonics of the 60 Hz line voltage. For reasons 
of stability, the modulation and reference signals were 
derived from.a {3 MHz) crystal oscillator. °The precise 
modulation frequency was 78.125 Hz, this being obtained 
Dyedigital, frequency «division of the: 3,MHz. signal. 

The input stage is common to both lock-in 
detectors. The secondary of the input transformer 
provides a balanced output which is connected to a 
differential-input amplifier. This arrangement pro- 
vides isolation and common mode rejection. The 78 Hz 
(w) Sinise” liste. Ve. (20) Signals are amplified and passed 
through a band pass filter, before being fed to the 
separate channels. 

In order to avoid overloading the later stages 
of the lock-in detectors with wide-band noise plus 


60 Hz or 120 Hz interference signals, each channel is 
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equipped with a band pass filter of bandwidth 10 Hz 
eentcredwabouts/csor 156°HzZ. A fairly close approxi- 
mation to a rectangular response in the frequency 
domain was obtained using multistage stagger-tuned 
active filters. Provision is made for the system to 
control the gain of the amplifiers, 

fnvour correlators, =the incoming Signal is 
multiplied by a square wave at the reference frequency 
Of 78 or 156 Hz. This is achieved by means of a 
digitally programmable amplifier the gain of which is 
plus one for the first half of the reference cycle and 
minus one for the second half. 

Two-pole Tchebyscheff low-pass filters, the cut- 
off frequency of which can be controlled by the system, 
remove high frequency noise components from the corre- 
lator Foutput:  thevoutput fof “the “filter is then! féd ‘to 
the analogue multiplexer. 

Diode shaping networks and band pass filters are 
used to convert the modulation signal to a sine wave. 
This is then applied to the modulation coils via an 
ac coupled current source. The ac coupling ensures 
that no de currents bias the steady magnetic field, while 
the use of a current source rather than a voltage source 
provides a more stable field modulation amplitude and 


phase. 





The phase of each reference signal relative to 
the modulation, is manually adjustable over a range of 
Poe eeeiteCanmaLSOsbee changed. bysexactly: 90° ‘and 180° 
by the system. 

The On and 20 correlators and the modulation 
source were constructed in a separate rf shielded box 
(called Superbox for reasons which will become clear 
later!), to prevent the mixing of the spectrometer and 
reference frequencies. This will be explained in more 


detail in, section, 5.4.3. 


5.2.12 Analogue multiplexer and A/D converter 


The various signal voltages to be measured are 
connected to the A/D converter via the analogue multi- 
plexer, which has fifteen differential input channels, 
Wren provision fOr gain auto-ranging Onsynine of them. 
Digitally-programmable amplifiers form the basis of the 
multiplexer, so that any one of the input channels may 
be connected to the A/D converter on command from the 
system. 

The A/D converter is an Analogic "three-and-one- 
half-digit" panel meter, model number 2532. The 
display gives a visual check on the operation of the 
equipment. For the A/D converter, a full-scale reading 


Of 22 ,000 counts) corresponds, tO an input of) +3 V.. The 
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analogue circuitry preceding the A/D converter can 
handle a voltage in excess of this, thus full use can 
be made of the converter's dynamic range. 

On completion of a conversion from analogue to 
digital representation, the count number plus the 
multiplexer gain can be viewed by the BAG) pe og 
computer via the monitor lines. The computer has a 
software subroutine for converting this output into 
the floating-point number format recognized by the 


computer software. 


5.2.13 Magnet and NMR field-lock 


The magnetic field Ho of approximately 6700 G 
is produced by a box-section permanent magnet with 8" 
diameter pole caps and a 2" gap. The magnet is 
"overshimmed" such that there is a field minimum at 
the center of the pole gap and a maximum at a radius 
of 2". The best field homogeneity occurs in the latter 
region and is about 100 mG over al cm volume. To 
minimize vibration problems, the magnet is mounted on 
antivibration pads. Thermal stability is achieved by 
surrounding the magnet with 1" thick styrofoam. Tests 
show that the thermal time constant of the insulated 
magnet is approximately 34 hours. 

The temperature coefficient of Hy was found to 


be -1.2 G/°C. Since we required a field stability of 
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a few tens of milligauss, it was necessary to stabilize 
Ho: This was achieved ister: a eh NMR field-lock, 
which provided a current fed to a pair of field correc- 
tion coils wound around the pole caps. The ava 
resonance from an aqueous solution of LiCl was detected 
using a Rollin spectrometer (Rollin, 1946) employing 
25 Hz field modulation and lock-in detection at the 
modulation frequency (see Fig. 5.6 for block diagram 
of the field-lock unit). The modulation signal and the 
reference signal for the lock-in detector were obtained 
by frequency division of a 1 MHz crystal-controlled 
Signal. The lock-in detector is similar to the Wn 
detector discussed in section 5.2.11. The 25 Hz modula- 
tion frequency of the Rollin spectrometer was chosen so 
as not to be harmonically related to the 78.125 Hz 
modulation frequency used in the Robinson spectrometer. 
The rf source for the Rollin spectrometer was a 11.1 MHz 
crystal-controlled oscillator of stability 1 part in 
10/7°c, The 11.1 MHz frequency was chosen so that the 
field correction would be both small, and always of the 
same sign so that a unipolar power supply could be used. 
In practice, the field correction was about -10 G. 

The Rollin lock-in detector provides a suitable 
error signal in the region of the resonance center 
(cf. the discussion in section 5.14 of the OW and 20 


lock-in detectors). Whenever this error signal is 
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The hlock diagram of the operation of the 
NMR Field-lock used to stabilize the 


magnetic field of the permanent magnet. 
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larger than a predetermined level, an interrupt signal 
is sent to the computer. This can then change the 
current feeding the correction coils thereby reducing 
the error signal to zero. This current is controlled 
by a 12 bit D/A converter, the bit size corresponding 
to 1.83 mG. The noise in the spectrometer output 
corresponds to less than half a bit, hence Hy is 
stabilized to a few parts in fooe 

The lock-in detector output is also zero far 
from the resonance condition! This POSSi DI Iity 1s 
detected by monitoring the output of a filter, tuned 
to the second harmonic of the modulation frequency. 
An interrupt is sent to the computer if the Output is 
too small, thereby alerting the operator that the field 
BLapilaZacton unit 1ss0ut, of Lock. 

The control algorithm for the "Field—-lock" is 
discussed in detail in section 5.3.6. 

tTheserteid-lock”™ actually stabilizes the field 
at the LiCl sample, whereas what is really required 
is to stabi la ze@thelim elds atithe NaNO, sample. This 
was achieved by positioning each sample on the median 
plane of the magnet gap and at the same distance from 
the gap center. For such a symmetrical arrangement, 
the correction field is the same at both sites. The 
samples were actually placed four inches apart, thereby 


ensuring that both were in the region of maximum field 
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homogeneity. If the changes in the field of the magnet 
are identical at each site, fthen Hy 1s also stabilized 
at the NaNO, sample. This was checked by changing the 
room temperature by 2°C and monitoring the magnetic 
meld lat theslicl and NaNO, samples. 

With the four-inch sample Separation and the 
40 mG modulation amplitude used for the field-lock, 
it was found that the 25 Hz fringe field at the NaNO, 
Sample was negligible. Likewise, the 78 Hz field at 
the LiCl sample was found to be sand However, the 
fringe field of the 78 Hz modulation will induce small 
voltages in the field-lock correction coils and the 
25 HZ modulation coils. To prevent these voltages 
being mixed with the small rf Signal which is invariably 
radiated by the Robinson spectrometer, passive rejection 


filters tuned to 78 Hz were included in the supply leads 


tothese coils. 


5.2.14 “Manual operation console 


This allows the user to control the system 
without recourse to the computer, this feature being 
very,.useful for maintenance and Lestingunpurooses:. 

The console consists of the following items:- an 
address select switch, sixteen data Switches, a three- 
position "function" switch, three "strobe' buttons, sixteen 


monitor lamps, various monitor outputs, and the digital 
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panel meter. The exact effect of the "function" and 
"strobe" buttons depends upon the particular device 
being addressed. However, the general procedure of 
data entry is similar to that used to the frequency 
control D/A converter, and we shall use this device 
aS an example. 

The address of the D/A converter is selected 
by the address select switch; the monitor and data 
lines are now connected to the device, and the monitor 
lamps display the contents of the D/A data latches. 

To change their contents, the new data are fed into 
thessi xtecnsdatasswitches, the, four.most,. significant 
Switches being in this case redundant. The function 
switch is then set to LOAD, the relevant "strobe" button 
is pressed, and the data are transferred from the data 
switches to the D/A data latches. The monitor ‘Lawes: 
now display the new data. 

It is sometimes desirable to be able to incre- 
ment or decrement the number stored in the data latches. 
This can be achieved via the console by setting the 
“Eunction" switch to UP. Each time the "strobe" button 
is pressed, the number is incremented by one. Similarly, 
when the "function" switch is DOWN, the number is decre- 
mented by one. 

All the functions which can be accomplished by 


hand via the manual operation console can be performed 
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much more rapidly by the computer when it is Conero ling 
the system. "7 in*this* mode,” the computer disables the 
manual operation console. However, since the computer 
is sending data to the device for only brief periods, 


manual operation is still possible most of the time. 


Deae bo COMpULer, sinterftace and teletype 


the computer, which is the ALPHA 16 model, 
manufactured by Computer Automation Inc., is a sixteen- 
bit stored program machine with 4 K words of magnetic 
core memory. The following processor options were 
added: a buffered teletype (TTY) interface, a "power 
fail restart" and a crystal-controlled real time clock 
(RTC). The TTY interface connects the computer to the 
TTY, thereby enabling data transmission between the two. 
The contents of the core memory are not lost during a 
power failure and the "power fail restart" option 
provides an orderly shut-down and restarting procedure 
if such an event occurs. The RTC provides, with the 
necessary software, the capability to perform operations 
in real-time. 

The interface is an in-house constructed device 
which is physically located in the computer case, and 
which is connected to the system by two fifty-conductor 
cables. Its function is to enable data to pass between 


the computer and the spectrometer system, while 
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allowing each to perform at their appropriate speed. 

PiGe bY £.:Si, a ect ASR33> which consists of a 
keyboard, printer, and a paper tape punch and reader. 
The latter feature provides a comparatively quick and 
easy method of storing and entering programs. 

We shall now describe the software aspects of 
the computer. Whereas there are many data analysis 
languages such as FORTRAN, there did not appear to he 
a widely accepted high-level real-time language 


suitable for our purpose. Therefore, one had to be 


created which would fit into the 4 K of memory available 


The alternative procedure of writing all instructions 
in machine language would, in the long run, have been 
much more time consuming. | 

There are basically two parts to the software: 
the control program which occupies approximately 3 K 
of memory, and the user program which can occupy the 
remainder. The control program is made up of a number 
of segments which collectively control the spectrometer 
system and TTY via the user program. Also in this 
part is a set of instructions used for "debugging" 
programs. The user program consists of a series of 
user instructions plus space for data storage and 
messages. Available to the user program are two 
working registers, the math register (MR) in which 


mathematical operations are performed, and the input/ 
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OutHuUL regrster (1/0 Reg), in which logical operations 
are performed and through which pass spectrometer 
System data and instructions. User instructions are 
in general two sixteen-bit computer words long. The 
first word is the instruction-class code while the 
second often contains an address to find or place data. 
Also included in the second word are up to four bits 
Specifying Which instruction within the class is to be 
performed. The instruction class set is as follows: 
MOLCeGeOne tii, AGG, Subtract, Multiply, Divide, Modulus, 
Jump, Increment and skip if 0, AND I/0 Reg, Complement 
170 Reg, Compare 1/0 Reg-and skip if’equal, Set RTC, 
Pescenle, shead A/D,ePack A/D reading. (This Jast 
instruction converts the A/D reading into the floating 
point format, mentioned in section 5.2.12.) This group 
of instructions enables the spectrometer system to be 
Grrr crenctysconcro. lea. Additional instructions could 
be added, but at the expense of the already small space 
available for the user program. 

Fach user instruction actually requires many 
machine instructions in its execution, a task handled 
by the control program. The control program also 


performs numerous other housekeeping tasks which are 


94 


invisible to the user, eg keeping track of time, handling 


interrupts, and checking user instructions for validity 


This last task becomes very visible if any errors are 
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detected, since the user program is then terminated, 
and a message indicating the type of error is printed 
onrthesity! 

The use of subroutines facilitates the writing 
of complex programs, and often enables them to occupy 
less memory space. There is provision for four levels 
of subroutines in the user program, this number being 
the result of a compromise between user language flexi- 
bility and available user language space. 

Both integer and floating-point numbers are 
available to the user program. fnLeger numbers occupy 
one memory word and have a range +16,383. Floating- 
point numbers occupy two consecutive memory words and 


ae 
have a range #10734 


WitLieanmaccuiraeylotvsevenldecimal 
digits. Floating-point numbers are used whenever a 
large dynamic range is important; integer numbers are 
used in cases, such as manipulations to generate computer 
addresses, where round-off errors are intolerable. 

In our real-time system, the handling of inter- 
rupts from the system is a vitally important aspect of 
its operation. We shall now describe this in some 
detail. There are two interrupt lines to the computer, 
one dedicated to the A/D converter, the other to the 
remaining five system interrupts. This second line is 


fed from an interrupt receiver which manipulates the 


five interrupts. To avoid conflicts of interest if two 
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of these interrupts occur simultaneously, they are rated 
in priority so that the one with highest DLLOristy 1s 
serviced first. In order of decreasing DIuLOuLty aac ne 
device interrupts are frequency counter, misbalance 
detector, and those for the NMR field-lock. The inter- 
rupt receiver stores each Gevice*interrupt an a iatch, 
the output of which passes through a number of masks 
before being connected to the computer interrupt dine, 
The latches may be reset, and the state of the masks 
(ie mask or unmask) set by the user program Via an 
Input/Output instruction. This rather complicated 
arrangement is essential because the computer may be 
performing a task in which it would be very inconvenient, 
Or even disastrous, to be interrupted by a particular 
device or devices. When an interrupt is received by 
the computer, it determines from which device it came 
and then jumps to the relevant user-program interrupt 
subroutine. 

An example Of a Situation an -which interrupts 
could cause a problem, is in the control of the A/D 
converter. JAfter -the computer Sendsea signalbetoetne 
A/D converter to start its measurement, the COMpucEeT 
waits for a signal on the first interrupt line which 
Signifies the end of conversion. If, for example, 
the field—-lock interrupt is*®received during this wait, 


the computer would service the field-lock. However, 
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the field-lock interrupt subroutine involves the A/D 
converter! To avoid this problem, the control program 
automatically masks the second interrupt line during 
the wait, and resets it to its original state after 
completion of the A/D instruction. 

The other major real-time function is the RTC. 
The user program sets the RTC in units of 10 ms by 
the set RTC instruction. A test RTC instruction checks 
Pearle slimesiseup,, and, if 1t is, the computer jumps to 
a specified program step. When going from a lower to 
a higher level subroutine, time is still kept in the 
lower level. It is also possible to set and test the 
RTC in the higher level subroutine independently of the 
lower one. 

Finally, we give an example of a Simple program 
which could be used to produce a linear field sweep. 
In this example, the field is not locked, and the sweep 
is produced by incrementing the contents of the field- 
lock D/A converter once per second. The program is 
Shown in Fig. 5.7, the numbers being in hexadecimal 
format in which decimal numbers ten through fifteen are 
written A through F respectively. The control program 
executes step O and proceeds @withwstepo el, 2,5 ,etce 
unless directed otherwise. The starting location and 
the storage space indicated can actually be anywhere 


within the user area of the memory. The parameters of 
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A sample program written in the user 
language. This program could be used 
to produce a linear magnetic field sweep. 
The “contents. column wa Wie benen tered 


into the computer memory. 


MEMORY 
LOCATION 
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the sweep can be changed by simply feeding different 
numbers into locations OEO9, OFOO and OFO] of the 


computer memory. 


52.16 System monitors 


One way of checking the operation of the hardware 
is to use the manual operation console to view the 
state of each device. However, in order to view para- 
meters stored within the computer without interrupting 
the experiment, two monitoring devices are connected 
to the Buss. Each contains a sixteen-bit data latch 
and an associated D/A converter; the user program can 
output important parameters to either of these devices. 

Other monitors which are not connected to the 
Buss are a Hewlett-Packard (HP) 1201Bs0oscilloscope to 
view the audio output of the Robinson spectrometer, a 
Darcy .TSI 385R frequency counter with decimal display 
connected to the output of the frequency converter, 
and two HP 680M strip-chart recorders to monitor the 
outputs of the Orn and 204 lock-in detectors. The event 
markers of these recorders are under the control of 


the system. 
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dete t/n, POWer Supplies 


ASSWr le beesecenein Section 51490, atmajor 
problem was the elimination of systematic errors 
arising from the interaction between the rf and audio 
Signals. The ac power line, ground connections and 
common dc power supplies provide a very effective means 
of coupling the two signals. In order to reduce the 
problem, the units associated with the rf Circuitry. 
of the spectrometer were physically separated from the 
low frequency digital and analogue Circuitry. Moreover, 
the rf and audio sections were fed from "independent" 
115 V ac power lines. Extensive rf filtering was 
incorporated in the power line feeding the audio cir- 
CULCTY. | 

Grounding, as opposed to shielding, is only 
important for operator safety. Indiscriminate grounding 
of equipment can cause "ground loops" thereby creating 
troublesome spurious signals. The ground pins on the 
ac power plugs of the individual units were therefore 
not used. Instead a cold water pipe was used as a 
common ground point for the various parts of the system. 

To improve electrical isolation, the Robinson 
spectrometer has its own 200 V- 200 ma and -6.3 V power 
supplies which incorporated extensive rf and audio 
filtering. Since the system contains approximately 


forty cards of circuitry, it was impractical to provide 
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an independent dc power supply for each card. Instead, 
a few high-power supplies were used, and isolation was 
achieved by incorporating shunt regulators on all of 


the cards. 


5.3 User language computer programs 
5.3.1 General considerations 


The flow-charts for the User Programs are shown 
in Figs. 5.8 to 5.16. Although these are not written 
in the user language, they do display the Logical 
sequence of events of the actual computer program. 

The flow-chart is read as follows: start at step O 

and proceed with steps 1, 2, 3, etc., unless directed 
otherwise by a branch step. These are written in the 
form Of ajquestion. - If the answer isti:yes, then branch; 
if the answer is no, then proceed with the next step. 

Normally, the computer is running the "Data 
Coltection Program” and its subroutines, which are 
"NZ", "156" and "Resonance-Lock". When an interrupt 
is received from the system, the "Data Collection 
Program" is halted and the computer then jumps to the 
appropriate interrupt subroutine, namely "Frequency 
Counter", "Field-Lock”™ or, "Misbalance Detector". On 
completion of this subroutine, the computer returns 


to the "Data Collection Program" and continues where 
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Figure 5.8 


Flow cHart of ther “Datagcolieet ion 
program!. This program uses sine 
"N2", "156" sand St Resonance. lock. 
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The "Frequency Sweep Program" is used to generate 
a triangular frequency sweep for the spectrometer. It 
is used instead of the "Data Collection Program" when- 
ever tepisenecessany tosploteéor search fora -resonance: 

Of the total memory space of 4095 words, only 
1160 are available for the user programs. Of these, 
approximately 1000 are occupied by the 500 user ins- 
tructions and the remainder are usable for data storage. 

We discuss in the following sections 5.3.2 to 


5.3.6 the programs used for the enhancement experiment. 


Seo <20n These DatasCollection" ‘Program 


As discussed in section 5.1, the On lock-in 
detector maintains the spectrometer frequency at the 
resonance center and the 2u lock-in detector measures 
the resonance amplitude. The output of the 200 lock- 
in detector is shown for one "saturation eyciles min 
Pigs so. 2. ence Une start OfeeniSecyCle, cues OCk—141 
detector rererence has a phase 0° (say) and the satura— 
tion power is off. Data is collected for ten seconds 
and then the reference phase is changed to 180°, thus 
reversing the lock-in detector output. After waiting 
1.4 seconds to allow the output from Tchebyscheff low- 
Dass -lelter, tormcone to within O[0S¢ of ats ftinal= value, 


another twenty seconds of data is collected. tThe 
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reference phase is now reversed back to 0° and, after 
the 1.4 seconds wait, ten seconds of data are again 
collected. This completes one "phase-cycle", which 

is repeated two more times. The saturation power is 
then applied and the resonance amplitude thus increases 
rapidly. (The noise level also increases due to noise 
from the saturation source,) -A time of forty seconds 
is allowed for the spin system to reach equilibrium. 
Six "phase-cycles" are then performed. The saturation 
power is thereafter removed and the spin system relaxes 
back to the unsaturated state. (This recovery allows 
the value of Wy to be determined as discussed in 
section 4.2,) Forty seconds after removing the satur- 
ation power, three more "phase-cycles" are performed. 
The "saturation cycle" is completed by the teletype 
printing out the experimental data. Every fifth 
Becaburallonncycle | some diagnostic. information, is 
printed with the experimental data. 

The flow charts of the program which generated 
the sequence of events just described are shown in 
Figs. 5.8 to 5.10. This program also analyzes the 
data collected. In these flow charts, Nl is the number 
of times per quarter "phase-cycle" that the 20, Lock= 
in detector is measured with the A/D converter; N2 
is the number of times per quarter "Saturation cycle" 


that the "phase-cycle" is performed, and N3 is the 
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number of "Saturation cycles" performed between diag- 
nostic data printing. In the present experiment, they 
Stem UUM eat respectively.) | Ini Figs. 5.8 to 5.10, 
Dis the true peak resonance amplitude for one 
"phase-cycle". The quantity U, is the mean over the 
i'th "saturation cycle" of the D values obtained in 
the absence of the saturation power. Similarly Si is 
the mean of the D values obtained in the presence of 
the saturation power. 

We shall now describe some aspects of the data 
analysis. The output of the 20, lock-in detector in 
general includes a dec offset. This offset, however, 
remains constant on reversing the reference phase, 
whereas the resonance signal changes sign. Therefore, 
the difference between the lock-in detector signals 
obtained at the two reference phase settings is a 
measure of the true peak-to-peak resonance amplitude. 
Also, because equal time is spent collecting data at 
each of the two phase settings, and because the 
"phase-cycle" is such that the data for each setting 
is centered about the same point in time, D is unaffec- 
ted by linear drifts in the dc offset. Higher order 
adriftts do, however, 1nctrodiice-an- error in Dp, but a1t is 
unlikely that such drifts would occur consistently 
over many "phase-cycles". Therefore, the error intro- 


duced will statistically average to zero. The 
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"saturation cycle" similarly cancels completely any 
error due to a linear gain drift and reduces any error 
caused by non-linear drifts. 

Ideally, the phase’ and "saturation" cycles should 
be short compared with the characteristic time of any 
drifts. However, there are practical limitations as 
to how short these cycles can be. One is the forty 
second wait required for the spin system to reach 
Sour tibrrome che Shorcer they saturationacycie, , the 
greater the proportion of time spent in this wait 
period, therefore less time is spent collecting data. 
The other waiting period is the 1.4 seconds required 
by the low-pass filter to stabilize after a reference 
phase-change. This time can be reduced by increasing 
the filter bandwidth. This would mean that the A/D 
sampling rate would have to be increased proportionally, 
LOmavOldsalrasing problems (Carlson, 1975) and loss of 
data. Our A/D converter has a maximum practical sam- 
pling rate of 10 Hz, which implies that components 
with frequencies above 5 Hz must be rejected. This 
was achieved by setting the 3 db cut-off frequency of 
ther Tschebyscheff ‘filter*to ] Hz. Hencesformithe cycle 
shown in Fig. 5.2, data was noiieored for an acceptable 
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Although Ss. and U., the enhanced and unenhanced 
resonance amplitudes, have normally distributed errors, 
the ratio E; = S;/U; will in general have a skewed 
distribution with a long tail on the high-value side. 
We therefore take as our best estimate of the enhance- 
ment EAV = 2S; /2U; rather than E the mean of the indi- 
vidual E.- The difference between E and EAV depends 
upon the error in U;, and in our case was found to be 
approximately 1%. We obtain the random error in EAV 
from the scatter of E; about the mean, using the 


expression 


(SEAV) * = 3 { (3B) - (rE)? (5.5) 


B|F 


The quantity SEAV could alternatively have been calcu- 
lated from the errors in LU, and iS. . However, this 
has the disadvantage that one would have to take into 
account correlations in U; and Ss eaused: by, dain 
variations. The effect of such gain variations on 
the value of E; is, of course, minimized by the struc- 
ture of the "saturation cycle" as described previously. 
The values of Us, Sir Ey EAV and (SEAV) * are 


printed on the teletype after each saturation cycle. 


A two hour sample is shown in Fig. 5.17. 
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A two-hour sample of the data "printout" 
obtained during the steady state double 


resonance measurements, 


Us -§.611847E+868 S= 
EAV= +1.453256£F+85 


U= -5.367164E+968 S= 
EAV= +1.562361E+86 


Us -5 ABS IATE+2S S= 
EAV= +1.493165E+58 


Us -§.376785E+688 S= 
EAV= +1.491928E+88 


U= -5.584389E+868 S= 
EAV= +1.-481558E+66 
T= +66335 

#61391 

#61563 

+G5666 

+6566) 

+86684 


U= -5.384547E+998 S= 
EAV= +1.483662E+88 


U2 -5.A472248E+86 S= 
EAV= +1.,487498E+96 


U= -SeSi6847E+99 S= 
EAV= +1.49838045+66 


Us -5.341648E+686 S= 
EAV= +1.4995345+08 


Uz -§.396581E+88 S= 
EAV= +1.497321E+89 
T= +686666 — 

+61396 

+61596 

+86669 

+66002 

+66062 


-8.15S416E+66 
+2e3S4136E-87 


2 | 
u 


-8.3399316E+58 E= 
+1.261592E-63 


-7.974368E+868 E= 
+6-228553E-94 


-7.981L791E+80 E= 
+3 eS44892E-84* 


-$.867274E+69 E= 
+2.972830E-64* 


-$.826966E+838 E= 
+2 .881129E-84 


-8.2848495+96 E= 
+1.691312E-94 


-8.349989E+688 E= 
+1-386186E-84 


-8.391632E+668 E= 
+1.794556E-84 


-7-962885E+88 E= 
+1.588782E-84 


+1.453258E+08 


+1.553718E+68 


+1. 474498E+98 


+148 4513E+86 


+1-444612E+68 


+1.498741E+88 


+1.-513975E+68 


+1.513763E+68 


+1.575984E+66 


+1.4771385E+68 


Lupe) 
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5.3.3 The "Frequency Counter Interrupt" Subroutine 


The loops which control the Robinson spectro- 
meter frequency are shown in Fig. 5.18. The frequency 
counter is contained within the FAST LOOP, the function 
of which is to maintain the spectrometer frequency at 
a predetermined "Correct Frequency". The purpose of 
the SLOW LOOP will be discussed in section 5.3.4. 

The flow chart for the frequency counter inter- 
rupt subroutine is shown in Fig. 5.12. This program 
is entered when the frequency counter generates an 
interrupt signifying the end of the frequency count. 
Steps 9 andllof the program provide frequency markers 
on the chart recorder in the following manner. The 
"Correct Frequency" is stored in the computer as a 


11] one Ww 


binary number. Whenever the 64 Hz bit is one, a 
is fed to the chart recorder event marker. Correspond- 
ingly, when the 64 Hz bit is zero a "zero" is fed to 
the event marker. Thus, the event marker generates a 
Square wave of "period" 128 Hz, whenever the frequency 
sweep seteereetn BE being executed. Steps15 and 17 remove 
the ambiguity introduced by the 12-bit count number 
supplied by the frequency counter. If the apparent 
change in the spectrometer frequency is greater than 


200 Hz, it is assumed that the frequency counter has 


malfunctioned for this count; therefore no correction 

















veoh, acer Soret xy * sit | on Be 


ue « oN welt fakdoes @oktw sqookl edt wit 


= 


'),2 ($28 gb Avods S26 yonaupsz2 sedam 5 

1 etd? abtltiw So wlnrau @£ 2a2n9oR) 
rogonorjosge stg Bliegaiee of af dokde 38: 
a an 

eft ."veasepag?s Jcats69" ban. Lereranees 6. 
ssc nt Eoegure!A of Iilgw S008 woe aaa 
past edd stot Seed 12>. eff? ¥; 


5 r 
-2% at swadea el gal sverdya Jae 


~< 
- * F A 1. 
; jad yonoeped?’ athe: aay Sides ei 
- 
3 om © . 7 ba. cc) oe > ci Sir ; ' ~ ee > eb ce a: tqexTSI8 


S Z 
Te 
f ypatt st ra me 2pow eit 40 24 Ons e aqedB 4) 
. } ; ‘<a 
7) aghageos oxtafs sas 


4 Dpeote er 7 VYormarpets SL Mb 


| Pé ed’ zeveret’ . .edaun — 


oo 






Teunem dndve aAbteoe? “teri, ety ot het) 
-* esse eboaper a Shae alte os 

te856 on tue @asvelantd seed -eaem 2K '® wide 
wWonkiter: of sovenedw \48 852 “be 28q" "20 “erm se) 
| reins ba Duce nek aie paiod-gt experi in 
Sar De 












= 


128 


= 
> -. 


— 4 ~. aA . 
rj: ) 1 IVE Le 


a) 


7 | : 4 
PA) 
7 





Pigure 2) 40 


Block diagram showing how the frequency 
of the Robinson spectrometer may be 


Varied ands CcOnLrLO.led,. 
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is made. Steps 25 and 27 check if the D/A number 
exceeds the allowable range of the 12-bit frequency 
control D/A converter. Steps 31 and 33 are used to 
count the number of times this subroutine is used. 
These numbers are part of the diagnostic information 


printed every five "saturation cycles". 


59.3.4 The "Resonance-Lock" Program 


If the resonance-lock is achieved in the usual 
continuous analogue manner, a systematic error in the 
determination of enhancement could result for the 
following reason. The output of the 20, lock-in 
detector is an even function of the frequency deviation 
from the resonance center. Therefore, the mean detec- 
tor output will depend not only upon the resonance 
amplitude but also on the magnitude of these frequency 
deviations. As the signal-to-noise ratio of the error 
signal is, in general, different for the saturated and 
unsaturated parts of the cycle, the magnitude of the 
frequency deviations would also be different. Hence, 
the value of E obtained would depend upon the difference 
in magnitude of these Frequency fluctuations between 
the two parts of the saturation cycle. This system- 
atic error can be avoided by keeping the locking 
frequency fixed throughout one complete "saturation 


cycle", and only making corrections to this frequency 
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between cycles. 

The "Frequency Counter Interrupt" subroutine 
keeps the spectrometer frequency locked to the "Correct 
Frequency"; the "Resonance-lock" program ensure that this 
frequency coincides with the resonance center. The w 
lock-in detector provides the error signal for achiev- 
tnigeciliomcOngte von acm natcated in the SLOW LOOP of 
Fig. 5.18. The error signal is sampled every twenty 
SeCOnds maSmindtcatecaminwrig. 65.9. (ror convenience, 
this sampling occurs every time the 20 lock-in detector 
reference phase is reversed.) At the end of each 
Saturation cycle (i.e. once per ten minutes), the 
cumulative error (278) is used to determine the required 
change in the "Correct Frequency" by means of the 


relation: 
SCorrecclonmiin Hz)! = Kx78/U; (5 26) 


asvindicated in Fig. 5.11.. The parameter k is an ex- 
perimentally determined scale factor. The cumulative 
error (178) is divided by the resonance amplitude Us, 
so as to avoid having to redetermine k for each 
spectrometer observing power or crystal setting. The 
"Correction (in H2)"is usually not, more whan 20 Hz, 

so if a correction of greater than 50 Hz is indicated, 
it is assumed to be due to noise rather than a genuine 


shift of resonance frequency. If this situation occurs, 
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a correction of only 50 Hz is made and the user is 
notified by a star (*) appearing after the data print- 
out of this cycle. Stars occurring in many successive 


cycles usually indicate an equipment malfunction. 


a oan The "Balance Network" Subroutine 


This subroutine can be entered if an interrupt 
is generated by the misbalance detector, indicating a 
Significant NMR sensitivity change, as discussed in 
section 5.2.7. The purpose of this subroutine is to 
readjust the balance network so that the misbalance 
detector output falls to an acceptable level (corres- 
ponding to an NMR sensitivity change of no more than 
0.01%). This is the "OK" level mentioned in the flow- 
Clartagiven ink iq ao. 1 47s Untortunatcelypzuthe misbalance 
detector gives only the amplitude but not the phase of 
the misbalance. This deficiency can be overcome by 
the following procedure: We measure the misbalance 
amplitude Ay and sehen turn-on. £or a .given. time, the 
motor driving the in-phase (X) mode. This changes the 
misbalance by a previously calibrated amount Ax in the 
x direction. We now again measure the misbalance am- 


plitude A and the x-coordinate of the final misbalance 
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This procedure is shown in Fig. 5.19. Once the x- 
coordinate is known, the X-mode motor is turned on 
for the correct time and in the correct direction to 
make the x-coordinate of the misbalance equal to zero. 
A similar procedure can then be performed for the 
orthogonal (Y) mode. This procedure is repeated if 
necessary, until the desired misbalance level is 
reached. 

For the above procedure to work in practice, 
a number of imperfections have to be accounted for. 
The output of the misbalance detector is not perfectly 
linear with input amplitude; this output is therefore 
linearized by a simple polynomial transformation. 
Also, the inertia of the motors and the backlash in 
the reduction gear-train are not negligible. The first 
“problem is overcome by waiting three seconds after the 
motors are switched off before measuring the misbalance. 
For the standard Ax or Ay shifts, the backlash problem 
was avoided by choosing the sign of the shift to be 
such that the motor turns in the same direction as 
before. Whenever a correction calls for the motors 
to turn in the opposite direction to its last movement, 
the motor is activated for a sufficient time to take 
up the backlash. Since a perfect backlash correction 


is not possible, the misbalance amplitude is remeasured 
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Vector diagram illustrating the 
procedure used to rebalance the 


"Balance network". 
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and the new x (or y) coordinate is calculated. A 
final problem is the Joule heating of the delay line 
which causes the orthogonal mode to drift for several 
seconds after the saturating signal is applied. The 
misbalance detector interrupt is therefore masked for 
the first thirty seconds to allow thermal equilibrium 
to be achieved and to prevent needless balance-network 
Gorrections. 

The correction procedure works quite effectively 
in practice, and only one or two corrections are 
required to achieve balance. The overall Stability 
of the network is such that at the minimum separation 
of 20 kHz between center line and satellite, the sub- 


routine was used only about four times persnour. 


5.3.6 The "NMR Field-Lock Interrupt" Subroutine 
ae ees E : 


This subroutine performs the software tasks 
involved with maintaining a constant magnetic field. 
The block diagram of the NMR field-lock is shown in 
EaAg.99.0) and che tlow-chart as given in Figs 65.13. 

A field-lock interrupt “is-generated Gf the Rollin lock— 
in detector output is greater than the deadband, or if 
the second harmonic monitor output is too small. The 

deadband voltage was chosen to be equivalent to a one- 


bit change of the field-lock D/A converter. This 
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subroutine checks if the second harmonic signal is 
PooeStaliancteoes Or, tlowechart)-; if at i¢, it “outputs” 
the D/A number to the field-lock D/A converter. This 

is done in case a stray pulse has altered the D/A 
contents. In step 9, the computer decides which Rollin 
lock-in detector interrupt was generated. It then in- 
crements or decrements the D/A number as appropriate. 
Step 13 is part of the diagnostic information printed 
every five saturation cycles. 

Because the magnetic field drifts are typically 
very slow, the "one-bit-at-a-time correction algorithm" 
was Satisfactory. Occasionally, rapid field changes of 
up to fifteen bits occur when iron objects are moved 
within the magnet's fringe field. However, even these 
Cranerentserrors are within the required stability of 


a few tens of milligauss. 


5.4 Systematic errors in the determination of the 


enhancement 
5.4.1 Classification of systematic errors 


There are basically three types of systematic 
errors in the determination of the enhancement; those 
relating to NMR, those relating to instrumentation, 
and those relating to data aquisition and analysis. 


Methods of overcoming systematic errors in the last 
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category have already been described in section 5.3.2. 
They consist of data aquisition in "phase" and 
"saturation’ cycles to eliminate offset and gain drift 
errors, and taking EAV rather than E as the best value 
of the enhancement. In the following two subsections 
we Shall describe the other classes of systematic 


CE EOrSi. 


5.4.2 NMR- related systematic errors 


We can determine W/W, using any pair of the three 


23ya mg9) NaNo.. There are definite 


advantages, however, to observing a satellite and 


resonance lines of 


Saturating the center line. As shown in section 3.2, 
the separation between the satellite and center line 

Po eCOmEt CGteOLdcL, = Droportional to the Z component of 
Enewsetectric=tield gradient at the nuclear site, whereas 
the much smaller shift of the center line is a second 
order effect. Lattice imperfections will cause local 
deviations in the electric field gradient, and it was 
found experimentally by Kornfeld and Lemanov (1960) 

that such imperfections result in a small loss of 
intensity of a satellite line, but a negligible loss 

in the intensity of the center dine. Lteis theretore 
safer to saturate the center line rather than a satellite 


to ensure complete saturation of the spin system. 


mee aun os 





Failure to achieve complete saturation would obviously 
lead to low aeheheoren: telhies. 

If any spurious resonances from the materials 
in the probe occur at the Robinson spectrometer fre- 
quency, they will introduce an error in the deter- 
Mination of the enhancement. Sodium is not a cons- 
tituent of any material used enn thes body sob sthe 
probe. However, since oe and oon have similar 
resonance frequencies, all coils and wiring within 
the probe were made of silver rather than copper. 


Another possible source of error is the NE resonance 


from the probe body. One would expect the Za 
resonance to be weak since the modulation amplitude 


Ins. sa Doses ea. Cade 


is very much less than the 
Psepanticularlystrue forfour present,.method of detection, 
since the signal is proportional to the square of 

the modulation amplitude provided this is small.) 
Moreover, the aluminium is at least half an inch away 
Promecherspectrometer, coil... iIn,addition, the.ri_observ- 
ing power used in the experiment is very much less than 
the optimum power for observing the tliat resonance in 
aluminium metal. A search over the frequency region 
where the resonances from NaNO. Can. occur, revealed no 


7 eee, 27, 
spurious resonances. In addition, the Al resonance 


occurs at a frequency lower than the aR center line 


in NaNO... We therefore generally made our measurements 
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on the high frequency rather than the low frequency 
Satellite, thereby ensuring that there was no possi- 
bility of a systematic error arising from ena 

For the enhancement value to give an accurate 
measure of W5/Wy, the resonance line shape must be 
unaffected by the saturation field. When the tran- 
sition rate caused by the saturation field becomes 
comparable with the spin-flip transition rate, 
"motional" narrowing of the line will occur (Bloch, 
beso seoatmles and COotts, 81958)". This would lead ‘to 
enhancement values which are too large. Estimations 
of the transition rates indicate that this regime 
occurs at a saturation level orders of magnitude greater 
than that used in our work. Also, no evidence of line 
shapes changes on saturation was found. 

According to the theory of Redfield (1955), 
Provatorov (1961)and Goldburg (1960,’ 1961), it is not 
possible to completely saturate a resonance by a 
steady-state technique, if the saturation frequency 
is mis-set from the resonance center. As we require 
essentially complete saturation of the center line, 
an Sane was performed to investigate how close 
to the center of the resonance the saturation frequency 
must be. By performing an enhancement experiment in 
the normal way, but mis-setting the saturation frequency 


from the center of the line, the results of Fig. 5.20 
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The enhancement of the high frequency satellite 
as a function of the frequency offset of the 
Saturation signal from the middle of the center 
line. The field modulation amplitude corresponds 
to 1.4 kHz in virequency) Unuts selhessol1 dameline 


is the best fit to the data points. 
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were obtained. The extended region of constant enhance- 
ment can be understood (Hughes and Reed, 1971) by 
noting ‘thatethe fieldymodulation has the effect of 
sweeping the saturation field back and forth across 
the resonance. Experimentally it appears that essen- 
tially complete saturation occurs provided the modu- 
lation takes the saturation field across the resonance 
Center vat sometstagervortthe imodulation cycle. 
Nearpethescnystal torientation (04= 9545, .the 
separation between a satellite and the center line 
is small, and it is possible for the "wings" of these 
resonances to overlap. Some of the saturating power 
would then "leak" into the satellite resonance by 
Spin-exchange interactions (Bloembergen et al, 1959), 
thereby altering the populations of the various energy 
levels. This point was investigated by performing an 
enhancement experiment in which the saturation power 
was applied to the wing of the satellite resonance. 
When the separation between the satellite and the center 
line is 40 kHz and the saturation power is set 20 kHz 
below the high frequency satellite line, we obtained a 
satellite enhancement E = 1.0025 +0.0040. When the 
saturation frequency was 40 kHz higher than the high 
frequency satellite frequency, we obtained E = 
1.0004 +0.0027. These results indicate that "overlap" 


effects are unimportant in our crystal, since we would 
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otherwise expect an enhancement of less than unity. 

This is another way in which the saturation field 
can introduce a systematic error. Sodium nitrate is 
not a perfectly loss-free dielectric (Mariani et al, 
1967). Thus, the large saturation field will directly 
heat the sample, thereby introducing an asymmetry into 
the two parts of the "Saturation cycle". This will lead 
to values of enhancement which are too small. The 
measurements deseribed in the previous paragraph there- 
EOL eesnowmchat overlap elrects: and yheating effects are 
Trevor. cicant. 


The linear E-P. extrapolation is only valid 


bs 


for small values of Be as was found experimentally by 


bs 
Reed (1970). He observed deviations from linearity 
when the increase 6E in the enhancement caused by the 
finite Dae was approximately one tenth of the enhance- 
ment value. In our measurements, for safety, we res- 


tricted our oer range so that 6E/E was always 


Ss 

less than 0.05. The data obtained at 0 = 180° 

shown in Fig. 5.21 show no evidence of any non- 

linearity, although this angular setting is one at 

which the satellite line saturates comparatively easily. 
As will be seen from section 6.4, the observed 

Signal decays exponentially with a time constant of 


approximately six seconds after removing the saturation 


power. To avoid the systematic error associated with 
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Figure 5.21 


Enhancement of the high frequency satellite 
as a function of the spectrometer power pamper 
The data were taken when the center line was 


severely saturated and the crystal orientation 


corresponded to 9 = 180°. 
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recording the unsaturated resonance amplitude while 

it is recovering from saturation, we must wait several 
time constants before recording any data. This is 

the reason for the forty second wait in the "saturation 
cycles". Calculation shows that at the end of this 
period, the unenhanced amplitude is 0.06% too large, 
and the mean error over the whole unsaturated part of 


the cycle is only 0.0013. 


5.4.3 Instrumentation-related systematic errors 


One error already mentioned in section 5.2.7 is 
the NMR sensitivity change Paced by a spurious injec- 
tion of saturation voltage into the spectrometer. 
Associated with this sensitivity change is a pulling 
of the spectrometer frequency (Hughes et al, 1973). 

In fact, for large injections, the frequency pulling 
can be so severe that the spectrometer is locked to the 
saturation frequency! Without the frequency-lock, the 
free-running spectrometer frequency will in general be 
slightly different for the two parts of the"saturation 
cycle". Data would therefore be gathered from two 
different points on the resonance. However, the 
frequency-lock ensures that the same point on the 
resonance is observed throughout the complete "satura- 


tion cycle’. 
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Another problem is the mixing of the modulation 
and rf signals within the spectrometer to produce an 
rf signal which is amplitude modulated at W, OF 200 
This will then cause spurious outputs at the lock-in 
detectors, whose signs are reversed when the reference 
phased is reversed. The modulation at 20 Teepe hit Cu= 
larly troublesome since it contributes a systematic 
error to the measured NMR amplitudes. A typical modu- 
lation depth for the NMR signals in our experiment is 


TOe2 


; we regard a systematic error of 0.1% of this is 
troublesome. This translates into a maximum spurious 
modulation amplitude of only 100 picovolts across the 
sample coil. 

One method whereby an amplitude modulation is 
produced is via the magnetoresistance of the silver 
Spectrometer coil. | This has thevetfectsnof directly 
modulating the Q of the tuned circuit. Indeed the 
only difference between this magnetoresistance effect 
and NMR absorption is that the latter is resonant 
whereas the former is not. However, measurement” shows 
that the magnetoresistance of silver is small, and it 
is estimated that a two gauss modulation field would 


12 


produce a negligible modulation depth of TOF at 2W 0° 
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We are grateful to Dr. J. Rowlands for performing 


the magnetoresistance measurement. 
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An amplitude modulated signal can also be caused 
by the non-linear mixing of the rf voltage and any 
voltages at War OF indeed 20 which enter the spectro- 
meter. Voltages at Ww, are more likely to be present 
than those at Zu, because of the low harmonic content 
of the modulation signal. Initial tests without the 
NMR sample showed that the U; and Ss values were non- 
zero, thus indicating an amplitude modulation at 2W0 
These U; and Si values were) founds to be roughly pro- 
portional to the rf amplitude and to the square of the 
field modulation amplitude. This indicates that mixing 
was occurring via a cubic term in the transfer charac- 
teristics of the rf amplifier. We now discuss various 
ways by which the modulation voltage can enter the 
spectrometer, and methods used to:overcome this problem. 

The field modulation will obviously induce a 
voltage in the spectrometer sample coil. The coupling 
to the rf amplifiers was therefore minimized by insert- 
ing a high-pass L-C filter between the sample tuned 
circuit and the first rf stage. The capacitor was 
approximately 3000 pF; the inductor was self-resonant 
at 2MHzZ, ‘and had a resistancesoft 5 2. sithis falter 
behaves like a 10°:1 voltage divider at 78 Hz and like 
a low-loss capacitative divider at 7.5 MHz. 

The fringe field of the modulation coils will 


also induce voltages in the spectrometer wiring. This 
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field was reduced an order of magnitude by placing a 
sheet of 1" thick aluminium between the spectrometer 
and the modulation coils. Further reduction was pro- 
vided by placing soft iron tube caps over the sensi- 
tive rf amplifying tubes. 

The electrostatic coupling between the modulation 
coils and the probe can give rise to small currents at 
Wn flowing in the probe body! These were minimized by 
feeding the modulation coils from a source which was 
balanced relative to ground. 

The modulation voltage present in the input tuned 
circuit can also produce an amplitude modulation by a 
totally different mechanism. This voltage will alter 
the resonant frequency of the tuned circuit via the 
variable capacitance diode. A frequency dependent Gani, 
for example, will then convert this frequency modulation 
into an amplitude modulation. The modulation voltage 
across the diode was therefore reduced by using four 
diodes in parallel and decreasing the size of the series 
Capacitor to the smallest value which would give the 
desired tuning range. As seen from Fig. 5.5, these 
capacitors form a voltage divider, which in the final 
configuration provided an order of magnitude reduction 
of the modulation voltage pickup across the diodes. 

Precautions were also taken to prevent signals 
at the modulation frequency from entering the spectro- 


meter by any other route. All connections, in 


ES : ” 
« 7 : $ ‘ 
s L. » ath é 
= 
7 \q 3 j 
‘ a 7 > lie 
oe . 
4 ¢ 
= k Z ss - 
‘pe a2 « 
‘ = 
. - 
a 
"ee = ae ’ f ewe 
’ s) | ae 
-* - 
- q 
i { uJ! ed 
i 
‘4 A . = 
-t 
: 
i | 
‘ 
~ 
— i 
7 ay™ - 
“ = ae 
i 
- a * ‘ 
i peur 
> ee 2m 
‘ .& La) 
. 
= : e, ' , 
f= bs < Pa 
wat s : , 
‘ “ee o — . a 
- aad > 
- . ~ ee 7 - 


; -> ~ 
. S 7 i 1G ? fi eae Os FY » = 
t 7 A se | 
—_ = rv, — ; 
a ia eS a 7 





=) ‘ te > 7 2 
. 2 wtf ss of —~" ; 
‘= - r. | 4 a y 4 i Vv oe == Ls i bs eA OF: 2106. a! : 
Th 
- ¥ a 








—— i. 
; = A ss. 


ital 


particular those from the power supplies,were exten- 
Sively filtered. 

Disappointingly, in spite of all these precau- 
tions a significant modulation was present even with 
the modulation source disconnected! This indicated 
that the problem was more subtle than appeared at 
first sight. It was finally traced to the following 
mechanism: radiation from the spectrometer mixed with 
Signals coherent with the 20, reference voltage and 
the modulated rf signal was reradiated back into the 
spectrometer! To reduce this strange phenomenon, the 
radiation from the spectrometer was reduced by includ- 
ing rf filters in all leads connected to the spectro- 
meter, except the one to the frequency converter. (All 
Ofte these precautions also reduce «the sensitivity of the 
spectrometer to ease radiation.) However, the line 
from the spectrometer to the frequency converter still 
provided an unfiltered path into the spectrometer, 

This line was therefore fed by an amplifier with a low 
reverse gain, thereby reducing the chance of radiation 
leaking in by this route. The frequency converter, 
itself, was mounted adjacent to the spectrometer, and 
was encased in an "rf tight" steel box. At this stage 
of the development, the wiring to the probe and balance 
network was changed from the original coaxial to tri- 


axial cable in order to improve the shielding of these 
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units. 

However, even these additional measures did not 
eliminate the spurious modulation completely! Tests 
showed that the amplitude and sign of the spurious 
U; and Si were strongly dependent upon the spectrometer 
frequency and level. Indeed, it was finally realized 
that anything that changed the rf radiation field in 
the laboratory also changed U; and Si. The nature of 
the mechanisms was finally confirmed when the lock-in 
detectors and their associated equipment were located 
in a room two hundred feet away from the laboratory. 

In this configuration the spurious U; and Ss values 

were found to be negligible. This was obviously a 
rather inconvenient permanent arrangement and another 
solution was sought! The cure was to keep the rf 
radiation field physically separated from the circuitry 
(ea On and 20, voltages. This was achieved by 
reconstructing the correlators for the lock-in detectors 
plus the associated reference and field modulation 
source in a separate shielded box, which we called 
"Superbox". All leads entering this unit were exten- 
Sively rf and audio filtered, and the box was constructed 
Erometh/ 3 steel. "RE gaskets" of aluminium foil were 
used to make the box as rf leak-proof as possible. 


Furthermore, the circuitry was constructed from low 
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power logic to reduce the currents flowing at the 
troublesome frequencies. 

A final test of two weeks duration found no 
evidence of any spurious output, and it is estimated 
that any residual modulation must be less than one 
part in iMate 

As mentioned in section 5.4.2, measurement of 
the enhancement of a satellite when the saturation 
power was "off resonance" gave an enhancement close 
to unity. These measurements were performed when the 
separation between the saturation and observing fre- 
quencies was 40 kHz or less. We also measured the 
"off resonance" enhancement when this separation was 
150 kHz. Rather surprisingly, the enhancement was 
found to be roughly 0.25% low. A careful study of the 
Pp dependence indicated that this systematic error 


obs 


decreases with decreasing ee Hence the error in 


bs° 


the determination of W./W should be small since they 


1 
are obtained from an extrapolation of the E-Pops data. 
We attribute this frequency-dependent systematic error 
to the following mechanism. A small amount of radia- 
tion at the observing frequency finds its way from the 
spectrometer into the rf power source Supplying the 
Saturation field. Since the balance network is only 


perfectly balanced at one frequency, there will bea 


net coupling between the saturation coils and the 
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spectrometer coil at frequencies other than the 
saturation frequency. Thus, the spurious signal 

at the spectrometer frequency will appear as an 
additional feedback path. This will then change 

the spectrometer level in a way which depends on 

its phase and amplitude. The amplitude will depend 
upon whether the rf power source is on or off. Thus, 
there will be a difference in the spectrometer level 
between the two parts of the "Saturation cycle", 
though a change of 0.25% would be tooywsmall) to, detect 
directly with our spectrometer level monitor. Since 
the NMR signal amplitude is proportional to the spec- 
trometer level, this mechanism will obviously lead to 
a systematic error in the measurement of enhancement. 
From an analysis of the characteristics of the balance 
network, we would expect that the error would be pro- 
portional to the difference between the saturation and 
observing frequencies, although this point has not been 
checked. 

Another possible source of error in the deter- 
Mination of enhancement is non-linearity of any of the 
signal-processing stages. Since a typical modulation 
depth for the NMR is only Oper it is extremely unlikely 
that a significant non-linearity would occur in the rf 
stages of the spectrometer. The signal plus noise 


levels in the audio sections of the spectrometer is 
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typically a few tens of millivolts. Inspection of the 
tube characteristics indicates that the non-linearity 
over such a range is very much less than 0.01%. The 
remainder of the signal processing was performed by 
Operational amplifiers, which, because of the very 
large amount of negative feedback, are linear to a 
high degree. Thus, the only significant non-linearity 
will be caused by any stage overloading. The gain and 
bandwidth of the various stages were optimized to make 
full use of the available dynamic range. Also, the 
Signal and noise levels were kept much smaller than 
the dynamic range. We note that in the stages just 
PrELOrs co, che correlators where midmndie is "narrow- 
band", the noise has a Rayleigh distribution (Cohn, 
1964). In the stages following the correlators, the 
noise should have a Gaussian amplitude distribution. 
The system, from the output of the spectrometer to the 
computer was nevertheless checked for linearity. No 
non-linearity was found to within the accuracy of the 
K Five and one half digit” voltmeter used to measure the 
input signal. 

Finally, we consider any systematic errors that 
can arise from thesA/D converter. {Its linearity, as 
specified by the manufacturer, is one half a bit. For 
typical voltage levels measured, this corresponds to 


approximately 0.05%. Another error can be introduced 
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by the following mechanism. If a voltage occurs between 
Count. numbers neand wt) ,~ the first number will result. 
Therefore, on average, the count number will be one 

half a count too low. However, since we adjusted the 
overall gain of the system so that the voltage feeding 
the A/D converter was in the upper part of its dynamic 
range, this last systematic error is estimated to be 
less than 0.03%. The performance of our A/D converter 
was checked and was found to give an error of 0.2 

rather than 0.5 count, which therefore implies a sys- 
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5.5 Reduction of the random error in. the enhancement 


measurements 
woe Prela modu bation considerations 


Since we are only interested in resonance ampli- 
tudes as opposed to lineshapes, we adjusted the modula- 
tion amplitude to give a maximum output from the 20 
lock-in detector when the spectrometer is locked to 
the center of the resonance. This optimization was 
performed at each crystal orientation because of the 
Significant change in the resonance linewidth with 
crystal orientation. 

We note that the output of the spectrometer will 
in general also contain components at 4u ow, aes 2nw- 
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regained by using a lock-in detector for each additional 
harmonic. However, the improvement in the error in E 
would be rather small. For example, if we recovered 

a 40, signal which was 25% of the amplitude of the 20, 
signal, the error would be reduced by only 3%. A more 
promising alternative is to use only the 20, LOCK im 
detector but choose the modulation waveform such that 
all of the other components are zero. A series of 

tests in which a different amount of a 30 component 
was added to the modulation signal indicated that a 
small increase in signal amplitude of roughly 10% was 
possible. However, this feature was not incorporated 
into the system because of the time required to construct 
and test a working system free from systematic errors. 
Another alternative would be to use a bisymmetrical 
square wave field modulation of large amplitude, since 
this would produce a 20, Square wave at the spectro- 
meter output. The 20, component would be 37% larger 
than in the previous case. However, in this scheme 
transient effects would have to be considered, and so 


this point was not pursued. 


5.5.2 Optimization of the :datarcollection 


To take into account the perturbing effect of 
the observing power, ae on the measured value of E, 


we use a linear extrapolation of the E-P ob graph to 
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Peete equals 7ero. 4» Since W/W, is derived from the 
extrapolated value of E, we should collect data in 
such a way as to make the random error in this value 
as small as possible 

As shown in Appendix III, the optimum data 
collection procedure consists of spending one ninth 
of the time at Bane equals the maximum allowable value, 
and the remainder of the time at a value of observing 
power a quarter of the maximum. If instead of taking 
just two data points we were to undertake six measure- 
ments at ae =a pe Opes) 60 4/0, 5/68and 6/6 0f 
the maximum, and spend time on each point such that 
each has the same weight, then the random error in 
the intercept would be 1.7 times greater than that 
obtained by the two point method, Thus; the total 
time required for the second method would be (ea 
times longer than that required for the first to-give 
equal error in the intercept. 

In practice we took a number of different data 
points. However, the Lone values were in two groups 
close to the optimum values described in Appendix III. 


This was done to provide checks on the consistency of 


the data. 
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5.6 Experimental details of the transient experiment 


As described in section 4.2 the magnitude of 
Wy can be determined by recording the recovery of a 
satellite line after removal of the saturation power 
from the center line. Such an experiment was performed 
at 8 = 180° using the computer-controlled spectrometer 
system. As before, the resonance amplitude was measured 
using field modulation at Wn =e phase sensitive detec- 
tion at 2u The “Field-Lock" and "Frequency-Lock" 
programs (see sections 5.3.6 and 5.3.3) were utilized 
in order to maintain the resonance condition. To record 
the recovery of the satellite without distortion, it 
Was necessary to use” a detector bandwidth of 5 Hz. 
The signal-to-noise ratio was then so poor that signal 
averaging was essential. In principle this could be 
performed using the computer. However, there was 
insufficient spare memory,so the signal averaging was 
achieved externally in a Fabritek Signal Average model 
number 1062. The data acquisition of this unit was 
under the direct command of the computer, and the 20, 
lock-in detector output was sampled every 80 milli- 
seconds for a total of 512 data points per sweep. The 
flow chart for the user program for this experiment is 
shown in Fig. 5.16. Approximately twelve hours of data 


were collected at each value of Laparay and the average 
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recovery was numerically printed by a Fabritek high 
speed printer, model number 201. 

The optimum "two point" method discussed in 
Appendix III was not used since the linear extrapo- 
lation of the decay constant had not been previously 
experimentally verified. Also, since the decay had 
not been experimentally shown to be a single exponen- 
tial, the complete decay was recorded far out into 


the baseline for approximately seven time constants. 


Foie Crystal) information 


Most of the data reported in this thesis were 
obtained using a sodium nitrate crystal provided by 
the Harshaw Chemical Company, Cleveland, U.S.A. The 
crystal was in the form of a cylinder 10 mm in 
diameter and 9 mm long and the three-fold axis was 
parallel to the end face of the cylinder. A spectro- 
scopic analysis supplied by the manufacturer showed 
the following impurities: Al, 3 ppm; Ca, 8 ppm; Mg, 
Siopme Gus seeleeDOM. ends oi is. lppm._ Woes traces ofthe 
paramagnetic impurities Fe, Co, Ni and Mn was found. 
The crystal was annealed by the manufacturer prior to 
shipment. 

A few measurements were taken on a second 


sodium nitrate crystal also obtained from Harshaw. 
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This crystal was 15 mm in diameter and 11 mm long and was 


of the same nominal purity. 


5.8 Determination of the crystal orientation 


It is desirable to be able to vary the polar and 
azimuthal angles 6 and ¢ independently of one another. 
The crystal was therefore installed in a nylon sample 
Norderp tiesuch a way that 1t could be rotated about two 
mutually perpendicular axes. The crystal itself was 
mounted in a collar’ which could be rotated about a 
horizontal axis relative to the sample holder. The 
collar was equipped with a graduated scale so that the 
crystal orientation could be set to within one degree. 
The sample holder could be rotated about a vertical 
axis in the NMR probe. It was equipped with a circular 
scale and vernier which could be set to an accuracy of 
about a quarter of a degree. The crystal was aligned 
so that the three-fold axis was along the axis of 
rotation of the collar. In this arrangement rotation 
about the vertical axis changed ® but not $9, while 
rotation of the collar altered ¢ but not @. 

The angle § was determined from the quadrupole 
splitting (see eq. (3.1) with n = 0), while the angle 
@ was determined by the following combination of X-ray 
and NMR techniques. A precession photograph taken 


with the X-ray beam directed along the three-fold 
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axis* allowed the angle 6 to be obtained relative to 
Eheeenyvatatmeyv indersrOuanmacCcuracy, Ofe41° 8 The 
centrosymmetric nature of the X-ray pattern (Buerger, 
1963) did not allow a distinction to be made between 

¢ and $+ 60°. However, this ambiguity was removed by. 
monitoring the NMR linewidth of the center line as a 
fUNCcELOnROLeOK eBymcomparion.with eq. (3.5), the 
theoretical expression for the second moment of the 
center line, the correct x axis could be selected. 
Moreover, careful linewidth measurement at 9 = 60° and 
120° (where the amplitude of the cos 36 term is greatest) 
enabled the position of the x axis to be independently 


determined to within one degree. 


m . 
We are grateful to Dr. M.J. Bennettof the Department 


of Chemistry for performing the X-ray work. 
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CHAPTER VI 
RESULTS 


6.1 Orientation dependence of the quadrupole Splbigting 
eee ee ee Se eet PO LOS PLE CLAY. 


The electric field gradient at a nuclear site 
is characterized by the quadrupole coupling constant 
e*oq/h and the asymmetry parameter n, quantities which 
were Introduced in section’ 2.3. “It *follows from eq. 
(3.1) that the first order quadrupole splitting between 
a satellite and the center line for I = 3/2 is given by 


V4, 75 = #(e%0q/4h) (3 See Gl Seas cos 26). 


(6.1) 


When this work was begun, it was believed (Reed, 
1970) that there was a possibility that the crystal 
structure of sodium nitrate did not correspond pre- 
cisely to that predicted from X-ray measurements. [In 
order to verify that the one nuclei did indeed lie on 
a three-fold symmetry axis, careful measurements of the 
orientation dependence of the quadrupole splitting were 
made. In order to detect a very small asymmetry para- 
meterivit vis "clearsfromeq.’-(6.)P)"that*the satellite 
splitting should be measured in the vicinity of 9 =90° 
for different values of ¢. The low frequency satellite 


line in sodium nitrate was therefore plotted on a chart 


recorder as the first derivative of the absorption 
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line. This was done every 0.5° near 9 = 90° for 
Olea oe se andecG.5 >.) ain order ito, obtain a precise 
measurement of the quadrupole splitting the "Field- 
lock" feature was operative and a linear sweep through 
the resonance was provided by the "Frequency sweep" 
facia yecescribed in section 5.32).  Thesresonance 
center frequency was taken to be that of the cross- 
Ever spoint Of tne finstederivative signal. To avoid 
the systematic error associated with the time delay 
of the low-pass filter on the output of the On lock- 
in detector, resonances were recorded with both the 
frequency increasing and frequency decreasing. The 
mean of the two cross-over frequencies was taken to be 
the resonance center frequency. Another systematic 
error can be caused by a drift of the sample tempera- 
ture since the quadrupole coupling constant has a tem- 
perature coefficient of -0.33 kHz/°C (Andrew et al, 
1962). The sample was fairly well thermally isolated 
from the laboratory, so that any temperature drifts 
would be slow. However, to minimize any residual 
systematic error, measurements were made alternately 
at the-two settings 6 = 9-3.5° and 865°, 

The results of these measurements are shown in 
Fig. 6.1. The solid curves represent the best fit of 
the theoretical form (eq. (6.1)) to the data. We 


find that the maximum satellite separation from the 
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Figure 16.1 


Measurements of the satellite splitting near 
8 = 90° “for ¢7= —3. 5 [eand cou =1c or ogee ic 
solid lines are the least squares fit to 


the data. 
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From this it follows that 


n= "(0162 1/8) e 10" 


from which it appears that the asymmetry parameter is 
indeed zero at the 22 Na Sveecse ine sOodmumenitrate: 

From these measurements we can also obtain the 
value of the quadrupole coupling constant. Indeed, 
a more reliable value will in general be obtained by 
making measurements near 8 = 90° rather than near 6= OSs, 
Since any misorientation of the symmetry axis out of 
the plane of rotation will give rise to an error at 
the latter orientation. An alternative method of 
measuring the satellite splitting at intervals through- 
Qube Cie tiull Us02 range of 6, sutfers from the addi- 
tional drawback that small errors in the angular 
setting have a large effect on the satellite separa- 
Clon near oe= 45° andeesce 

The quadrupole coupling constant derived from 


Our measurements in conjunction with eq. (6.1)-.is 
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This value is in agreement with the value 
B04) 2ekh2 o1vensoverounds (1950); athe value of 334.2+0.8 
kHz deduced from Andrew et al (1962), and the value of 
335.4 kHz deduced from D'Alessio and Scott (1971). The 
fact that our value is marginally larger is not 
Surprising since the other measurements suffer from 
the drawbacks mentioned above, which would give rise 
to a slight underestimate of the quadrupole coupling 


constant. 


6.2 Orientation dependence of the linewidth 


Firstly, the ¢ dependence of the NMR linewidth 
was measured to remove the ambiguity in the X-ray de- 
termination of the angle ¢. Also in view of the 
anomalous $¢ dependence of Wo/W, reported by Hughes and 
Reed. (1971), a careful measurement of the linewidth of 
the center line and a satellite line was performed. 
The measurements were made at 9 = 60° and 120°, where 
the amplitude of the > dependent term is largest. 
Assuming that the linewidth is proportional to the 
square root of the second moment, it follows from 
eqs. (3.4) and (3.5) that the width of the center line 


is approximately given by 


AF A SS Le SOR COS oo ja KHZ (6,02) 


at @ = 120°. The corresponding expression for the 
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width Of a Satellite line is 
AF , =wvnoo9 (49 +50, L89ncos 30) kHz (oee3) 


at 6 = 60°. 

The peak-to-peak linewidth was measured by 
recording the output of the wn lock-in detector when 
a linear frequency sweep through the resonance was 
performed. The field modulation amplitude was typi- 
cally one sixth of the linewidth. 

TMYesdata are shown in ©igs.6.2 and 6.3.° The 


solid line is a least squares fit to the function 
A + B cos 36 (6.4) 


where ¢ is the value deduced from the X-ray work. The 


least squares fit gave the following results 
AF, = (1.755 4 0.010) [1- (0.269 +0.006)cos 396] kHz 


AF, eae O44 Oost (0. 200)7 05007) cos <30]* kHz: . 


ae nuclei 


The dataw clearly, confirm that the 
lie on a three-fold symmetry axis, and that the X-ray 
determination of the angle $9 is correct. 


The agreement between the theoretical and ex- 


perimental linewidths is surprisingly good, considering 


the approximation (eq. 3.4) that linewidth = 3 yM.. 


The agreement for both center line and satellite line 
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Peak-to-peak linewidth of the center 
line resonance for @ = 120° asa 
function of the angle 6. The solid 
line is a least squares fit to the 
expression 
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Birogunes. 6.55 


Peak-to-peak linewidth of the high 
frequency satellite resonance for 

8 = 60° as'a function of the angle 

?. The solid line is a least squares 
fit to the expression 
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indicates that contributions to the linewidth from 


lattice defects are negligible in our crystal. 


6.3 Orientation dependence of W/W, 





The E-Pobs data for each crystal orientation were 
fitted to a straight line by a weighted least squares 
fit technique. The error for each point was given by 
the scatter of the Es as discussed in section 5.3.2. 
The goodness of fit was not used to determine the error 
in the intercept because, for the small number of data 
potntsrused; | thisevalue is unreliable. ) Instead, the 
errorawasr calculated) using eq. (A3.3), an@ this error 
value will be quoted in all cases. 

The correction for the non-infinite saturation 
power was determined from measurements at each crystal 
Settingsusing equation §(4.8) or (4.9) asPappropriate. 
The value of the extrapolated enhancement obtained from 
the E-P obs data was multiplied in each case by the small 
correction factor € defined by 


E at infinite saturation power (6.5) 


E at the power used for the Pane data 


Cs 


The parameter 6, which is required in the calculation 
of Wo/W, was determined in each case from the quadrupole 


Soin cing. 
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The results obtained at § = 180° and db = -3.5° 
for Cases 1, il and Iii are given in Table L. The 
quantity E is the enhancement obtained by extrapola- 
ELONEEO Peete and Paneer The values of DS seace 


were used in conjunction with eqs. (4.4), (4.5) and 


(4.6) giving the following values: 


Bo 24 e270 O06 


= 
No 
ps 
= 
fb 
Il 


(6.6) 
O..0034240 20031 . 


= 
WwW 
Ss 
= 
a 
II 


The quantity W3/W) represents a small amount of mag- 
neticerc!laxation, occurring in the crystal, presumably 
due to paramagnetic impurities. 

It was shown by Niemela (1967) that W3 fox AS ee 
in sodium nitrate depended only weakly on crystal 
Orientation. Since Wi also does not: vary by more than 
415%, we make the simplifying assumption that y or | 
W/W, is independent of crystal orientation. A small 
Magnetic relaxation correction was therefore applied 


to all subsequent Case I measurements by using the 


relation: 
x= (E-1) /[2-E-6-(3y/(1+x))] . at G 7) 


Using the numerical values of Mi 333/111" 
M3333/@)411 and My433/"14111 given by Reed (1970) the 


derivatives 8 (Wo/Wy) /9 (My 343/My 471) etc. were 
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Table I 


Enhancement measurements obtained at 9=180° 


Case I 
E = 1.4894+ 0.0023 
6 = Pf70023 
ara tia 1492 8F 50). 0023 
Case If BS 1.504 0.0035 
b=) 0023 
Radars ee Peo hs62.0. 0035 
Case 7 lr. 
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Calculated and are shown as a Function Of£90 in Pig. 
6.4. These curves were used in choosing the 6 values 
at which the enhancement should be measured. For 
example, points near 6 = 60° or 120° are particularly 
useful in determining My43/q111° 

The results for the 6 dependence of Wo/Wy LOY 
Of= =S-5°eareushown vin Table I7T. Unfortunately, it 
was discovered that the trigger level for the Mis- 
balance Detector had not been set low enough when the 
Separation between center line and satellite Wasa 20m KHZ» 
When a separate experiment was performed to measure the 
NMR sensitivity change which had occurred under these 
circumstances, it was found that all enhancement values 
taken at a separation of 20 kHz should be multiplied by 
the small correction factor 1.0023. Measurements at 
- 40 kHz separation indicated that no Significant system- 
atic error was present. Each determination of W,/W, 
required approximately forty-eight to sixty hours of 
data, depending upon the Signal-to-noise ratio at the 
particular orientation. The results of Table II are 
Shown graphically in Fig. 6.5, the solid line being the 
non-linear least squares fit to the theoretical expres- 
Sion eq. (3.15). In order to estimate the goodness of 
fit we daeur ated the values of 2 per degree of free- 


dom,and the mean deviation defined as 


(6.8) 
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Figure 6.4 


The sensitivity Of W/W, tOmanVeLlatiton 


1313¢ M3333 or Mi413 as a function 


of M 


of the angle 8. 
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Table II 


The @ dependence of W5/W, at. i= —=o..5° itor Case: 1 


§° 


103.4 
sie ay. 
120.1 
129.8 
etal 
141.1 
148.3 
163.3 


WU 


E 


1.4893+0.0030 
1.497440.0033 
1.5069+0.0031 
1.5012+0.0033 
1.4742+0.0035* 
1.4695+0.0031 
1.4542+0.0028 
1.480640.0050 
1.5100+0.0022 
1.5306+0.0036* 
1.5367+0.0038* 
1.5416+0.0035 
1.5335+0.0034 
1.5208+0.0039 
1.5028+0.0040 


Mean of Cases I, 


lee 


POO LG 
ee Eo 
OOO 
0020 
1.0023 
O01 
T0017 
0016 
re OFLG 
i O06 
10016 
1.0025 
UO 
1.0024 
10020 


eee: 


a, 
Meooo oso 1 
1.052+#0.013 
1.082+0.013 
1.040+0.013 
0.928+0.012 
0.846+0.011 
0.867+0.010 
0.961+0.018 
1.075+#0.009 
ToT aoeue 0 te 
1.189+0.018 
122805017 
1.190+0.016 
eso “on cda 
TO Te CMON 


1.024+0.006 


- 
A "Misbalance" correction of 1.0023 has been applied, 
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Figure 6.5 


The orientation dependence of W/W, asa 
function of the angle 6 for ¢ =--3.5°. 

The solid line is a least squares fit of 
the data to the theoretical form. The 
dotted line represents the results obtained 


by Reed (1970). 
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where P is the number of degrees of freedom, and Xin 
PoMCnemhactwtUcuve MiclOtux.  hOrethe data. given in 


Table II we found 


me per degree of freedom = 0.49 


and Mean deviation OedijrsO0e28o. 


These values indicate that the theoretical form is an 
excellent fitting function to the data. The low value 


of “7 per degree of freedom (we would expect a value 
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close to unity) could be interpreted as either due to good 


fortune or to pessimistic error limits. One possible 


explanation is the fact that PS was found to vary up 


bs 
to +5% during the longer runs. Because of the non- 

zero slope of the B-Po, graph this will introduce an 
additional random error in the Ey values. This leads 
in turn to a pessimistic estimate of the error in the 


enhancement, since in practice, we recorded P at 


obs 


frequent intervals and used the mean value of Bene in 
subsequent analysis. In other words, eq. (5.5) gives 
a reliable error estimate only if the observing power 
is essentially constant. An estimate of the magnitude 
of this effect shows that it is capable of explaining 
the low value of x7 per degree of freedom. 

From the least squares fit to the data in 


Fig. 6.5 we found the following values for the fitted 


parameters 


“Oe 
* -% 










~ 2 a ‘; . 
i ’ 7 - 
‘ 4 ; a 
‘ > ' : iv vA 
be 4 i] 7 - i] 
ae 
; ‘ - me @ * =o ho = 7 ieee ort - 
<3 16 eesxpeb fo" sedime aaa eke Se 
' ; . n 
i _— “a == 1 Ms ics fp ir 
’ - - , Te ‘ 4 7 
. jee Ramey t a Fer we ~ _- oe id le _ 
; er 2 — 
> dy  * ee 
¥ ry 
> ; cae 
hws = i bG SsaiPplw ’ at 4 
= 
— + i “~ hE! <a x ‘e te * 
ba*. - Bf 3. 5a 
i 
a 2 = a “* -s~ = ~~ + a a = —— 
é : 4 4, -_— 
: | 3 “i |} 4803 -S236n L051 890.0" 
+ GF noisnast? earsss sas. 
“c)* de. 25 Tia oe ¢ Ye 
= 3 * > ~~ 2 ‘ 3 : =& £ ey Si t6q 
: * 
; % - Fé > _ bes 1G ee x =e ats Fi od | = am 
a bani he i nrg 
> ry smif worse tideinilegeaq-or, $e wIt 
a = ’ : 
‘ a0-.02-f tee $:d08t afeceh- nOtge 





eid. pas ~~ab- ee -, 


184 


My 3733/1111 =—0 20025 °%,0,0025 
M3333/Mj4q11 eee LOo 0 O00 59 (om) 
My 413/"1111 =—) L001 = .0.0022 


In section 7.1 we estimate the magnitude of systematic 
errors in the determination of the enhancement in 
order to make a realistic estimate of the error in 
the above quantities. 

It can be seen from Fig. 6.5 that there is a 
Significant discrepancy between the present results 
and those of Reed (1970). As mention in section 5.1, 
the method of Reed consisted of measuring the peak-to- 
peak amplitude of the absorption derivative when a 
linear frequency sweep through the resonance was 
performed. To check that the discrepancy between our 
results and those of Reed was not due to the difference 
in detection methods per se, we performed a "Reed type" 
experiment with the present apparatus at the crystal 
orientation @ = 180° where the discrepancy is largest. 
Repetitive linear frequency sweeps were performed under the 
control of the computer,and the output of the Wr lock- 
in detector was recorded in the Fabritek signal averager. 
The averaged resonances were then displayed on an X-Y 
recorder and the peak-to-peak resonance amplitude was 
measured. The data for this experiment are shown in 
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obtained using the "first derivative" method. 
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144.0 e500 005 
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and fee O015 
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in excellent agreement with the value P4928 450-0023 
found previously (see Table I). Thus, we conclude 
that the discrepancy is not due to the use of a 
different measuring technique. Other explanations 
for the discrepancy will be discussed in section thes oe 
Because of the anomalous 9 dependence of Wo/W, 
reported by Reed (1970), we undertook a detailed study 
Of Wo/W, as a function of ® near 6 = 60° and 120° where 
the amplitude of the $9 dependence is Tfargest.” *The 
results are shown in Table IV. Since the separation 
between the center line and the satellite line was 
20 kHz for these measurements, the previously described 
correction of 1.0023 has been applied. The corrected 
results are shown graphically in Fig. 6.6. The solid 
line is the theoretical form calculated using the M- 
values found from the 6 dependence. In other words, 
it is a "zero parameter" fit and there is no Wiehe) Cope 
any anomaly in the 9 dependence. The value of ye per 
degree of freedom is 0.73 and the mean CE tLOiy 
Ger iedet need.m( Gc) pss 0. 0.40110 29 We again find 
that Ee per degree of freedom is substantially less 
than unity, thereby confirming our SUSpLCion that the 
error limits are pessimistic on account of the varia- 
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Table IV 


The » dependence of W/W, aU oe 20 the srOreCcace al 


ed 
je9) 


* 


o° E gE Wo/W, 
eo) hee tO S0036 ie 0.U 6 1 MUSEO) ONG) 
PAs pS P4905 02020033 1.0018 ye Oe OO 01193 
Oreo 4d 2e2 0.0088 1.0026 OFS 0 On 
SG ao 1.4909+0.0034 1 SOMGREY iO OL 0.20 13 
Wea Gree eo 25520. 0030 OS Pal 4320 2013 
146.5 Mao O LE O0s0024 peu 0 le 0 452 O00 
dr, 1.4718+0.0039 0.0 29 OBIS 02040 be 
20625 1.4902+0.0036 PacOLs 1.000+0.014 
236.5, aoe. oe 020034 PaO0L9 ieee Or 0 15 
266.5 50082070087 0 0:E3 O42 Oe 
ea Gis 1.4660+0.0034 1.0028 ORM se O25 One? 
OO Bet 2 TSO 20035 Pe O019 ae O09 20m 4 


* 
A "Misbalance" correction of 1.0023 has been applied. 
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Figure 6.6 


The orientation dependence of W,/W, asa 
function of the angle » for 9 = 120.1°. 
The solid line represents the theoretical 


prediction. 
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the theoretdcal 


In order to check the consistency of our PaSilis, 
we have undertaken some measurements on the second 
crystal of sodium nitrate discussed in section Droid ts 
Results obtained at the orientation 6 = 0° for Cases 
I to IV are given in Table v. Usingregs. (4.4) ‘to 


(4.7), we obtain the values: 


W/W, = 1.020 +0.004 
(675710) 
= O.. 203 
W/W, 009750500 10 
Comparison of eqs. (6.6) and (6.10) shows good agreement 


between the Wo/W, results obtained with both crystals. 
The magnetic relaxation contribution W/W, is approxi- 
mately three times greater in the second Ccrystadl) than 

am the first. This presumably is an indication that 

the paramagnetic impurity concentration is not identical 


in the two crystals. 


6.4 Measurement of W, by the transient method 
ere eee en ee eee ee 


In order to obtain a value of Wi we measured 
the time dependence of the amplitude of the low fre- 
quency satellite at 6 = 180° after removal of the 
Saturation power from the center line (as described 
in section 5.6). This was repeated at various values 


of P - In orders that the data could be reasonably 


obs 


fitted to a single exponential (see section Aio2)) ale 
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Table: V 


Enhancement measurements obtained at 9 = 0° using the 


second sodium nitrate crystal. 


Case I 
E = 1.4859+0.0030 
ES=" EL0017 
E = 1.4884+0.0030 
corr 
Case II 
Be= 25015 26:20028 
So Jieaendy, 
E = 15504040. 0028 
CORY 
Case III 
E = 176568+0.0025 
cS Gk ALO Om Ry, 
E = 1.6596+0.0025 
Core 
Case IV 


E = 1.636740.0025 
Evel 0016 


E 2 Ps69S R50. 0025 
COE 
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the a ae values were small. As expected, the signal- 
to-noise ratios were poor, thereby requiring many hours 
Of data collection. A typical example of a twelve- 
hour run is shown in analogue- form iniFigq. 6.7. 

The usual method for analyzing an exponential 
Gecay as to take the logarithm of the data before 
Petting tO a) straight line.) The vale of the baseline 
of the data is then adjusted to minimize con This 
method is only applicable if the noise is small com- 
pared with the signal. A systematic error in the 
determination of the decay constant would otherwise 
Rect Since the error bars would be asymmetric after 
taking the logarithm. A more Satisfactory method is 
to use a non-linear least squares technigue by fitting 


the original data to the theoretical form 
Bete Beexp (=A tye 


In this method, there is no restriction on the signal- 
LO=norse ratio, of theydata, and no systematic error 
Should occur. Consequently, all the data collected 
is usable, including that in the "tail" of the decay. 
We used such a non-linear least squares fit 
method to analyze all of our numerical data, and the 
decay constant X is shown as a function of Pee in 
Table VI andi in) graphacas POrm in hig. 96.8", According 


to eq. (4.11), X should be a linear function of et 
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Figure 6.7 


The recovery to equilibrium of the high 
frequency satellite after removal of the 
saturation power from the center line. 
The trace represents 12 hours of data 


collection at a low value of P f 
obs 
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Table VI 


The decay constant A as a function of tes obtained 


at 06 = 180°. 
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and the solid line in Fig. 6.8 is a weighted least 
Squares fit to that theoretical form. The intercept. 
is equal to 


OeT622 08001] meecue 


In this case, the standard error is given by the 
goodness of fit, since the fit is not quite as good 
as one would expect from the errors of the individual 
data points. This discrepancy may be due to the 
approximation of fitting the decay to only one 
exponential, Or it may be due to small variations 

in spectrometer level during an individual run. 

The value of the intercept is used in conjunc- 
tion with eq. (4.12) to determine Wy: Using the 
previously determined value y = 0.0034+0.0011, we 
find 


(8.032+0.0084) x iM eee 


= 
I 


(65.41.15) 


4 cok 


W (2,740.9) x10 ~ sec ; 


3 


frspecte1on Orlcd.wio..2)eishnows that onlyeat) oa=—0sor 
180° is Wy a function of just one M-component. Thus, 


we can deduce immediately that: 


i sol 


262 = (4,819402050) 105" secre” . Comr) 


SN 314 


The smallness of the magnetic relaxation con- 


tribution can be appreciated by noting that if this 
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Figure 6.8 


The variation of the decay constant i as 
a function or Spectrometer power vena 
The solid line is a least squares fit to 
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were the only relaxation process, the spin-lattice 
relaxation time for the satellite line would be 


approximately 20 minutes instead of 6 seconds! 
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CHAPTER VoL 


DISCUSSION OF THE RELAXATION MEASUREMENTS 


7.1 Estimation of systematic errors in the relaxation 
data 

Rromelables Ji ands iVewe see that a typical random 
error in the extrapolated enhancement is 0.22. Thus, 
we need only consider systematic errors greater than 
about 0.01%. Of the systematic errors considered in 
section 5.4, we find that the only significant ones are 
of instrumental origin, which we shall now discuss. 

The trigger level for the "Misbalance detector" 
was set at a level which corresponds to no more than a 
0.05% NMR sensitivity change. (The only exception was 
when the saturation and observing frequencies differed 
by 20 kHz; a separate correction was applieds ior this 
case, as indicated in section 6.3.) We estimate that 
the systematic error in the enhancement caused by this 
source is about 0.03%, since the misbalance Signal is 
near the 0.05% trigger level only a small fraction of 
the time. 

Since noise is random, there is a finite 
pepe te of it exceeding the dynamic range of the 
signal processing circuitry, however small the standard 


deviation of the noise. As mentioned in section SLAY) Gane 
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this problem is most serious in the CrLCuner ver ust 
prior to the correlator in the lock-in detectors. 
Using the graphs given by Cohn (1964), in conjunction 
with a knowledge of typical Signal plus noise levels 
in our amplifiers, we estimate that the systematic 
error from this cause is no more than 0.052. 

The A/D converter in our system has a maximum 
non-linearity of one half a bit and a resolution of 
twelve bits. For the typical signal levels measured, 
this corresponds to a maximum Systematic error of 
approximately 0.058%. 

We now consider the magnitude of any residual 
Systematic error caused by the mixing of the rf and 
modulation related signals. The two-week check 
mentioned in section 5.4.3 revealed no Sionafticant 
Systematic error. Since a typical enhancement measure- 
ment required forty-eight hours of data, we estimate 
any residual error to be less than 0.35 (i.e. Y2/14) 
of the random error in an enhancement measurement. 
This corresponds to a 0.07% error. 

We have discussed in section 5.4.3 the frequency- 
dependent systematic error which we believe is due to 
stray signals at the spectrometer frequency being 
amplified by the rf power source and reentering the 
spectrometer via the saturation coils. Since this 


systematic error was observed to decrease with decreasing 
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Bowen the error in the extrapolated enhancement should 
be no larger than about 0.073%. 

The first two systematic errors discussed here 
will lead to enhancement values which are too low. 
However, the remaining errors could be of either S.0n. 
We therefore combine the errors in the usual way to 
find a probable systematic error in the extrapolated 
enhancement of about 0.15%. If this systematic error 
is uncorrelated from point to point, then the "goodness 
Ofeuit Weerror simi ts onaethe ratios of .the M-components 
‘Should include the effects of these uncorrelated 
Systematic errors. However, the systematic error Ih 
by its very nature, likely to be correlated from DOM te 
to point. To take this into account, we increase the 
error limits given in eq. (6.9) and our final estimate 


of the ratios of the M-components is 


My 31371117 = 0-8025 4 0.0035 
Me 353/My5 11 ie 1091 4.0.0083 (7B sI) 
Mi 313/83117 = 70-1001 20.0031 : 


The good agreement between the W/W values 
obtained with the first and second sodium nitrate 
crystals, indicates firstly that this quantity does 
not appear to be sample dependent. Secondly, since 


the second crystal is considerably larger than the 
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first, it shows that any systematic errors associated 
with sample size are negligible. Possible size depen- 
dent errors could be caused by large inhomogeneities 

in the rf fields or weak background NMR Signals emenat- 
Lnge from thes probe: 

We shall now consider possible systematic errors in 
the determination of the decay constant by the tran- 
Sient method. Unlike pulsed NMR experiments in which 
the spectrometer has to recover from an extremely large 
Overloading of the input stage, our method does not 
suffer from such-transients. Thus, shifts in the base- 
line of the data are not likely to be a significant 
problem. Also, the NMR information appears as a small 
modulation of an rf carrier. Thus, linearity problems 
in the detectors will not be significant. The finite 
bandwidth of the detection system can introduce a 
weeranat error in the determination of the decay 
constant. The low-pass filter following the 20 lock- 
in detector had a time constant of 0.03 seconds 
(corresponding to a 3 db frequency of 5 Hz). Since 
this time constant is so much shorter than the 6 second 
time constant of the observed decay, it is clear that 


this systematic error should be negligible. 
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Tee COomparisonmot the W5/W, results with theoretical 





predictions 


The excellent agreement between the 6 and ¢ 
dependence of W/W, and eq. (3.15) confirms the validity 
of the theoretical treatment of the orientation depen- 
dence of Wy and Wo. It also shows that the M-components 


which occur in Wy are equal to those in W Lnetusn, 


2° 
this implies that the spin-lattice relaxation is pre- 
dominantly caused by a multiphonon interaction. 

It was found (Satoh, 1965) that spin-lattice 
relaxation caused by vacancy diffusion becomes signi- 
ficant at elevated temperatures (>400 K) in sodium 
chloride. They therefore consider the possibility 
that. such a mechanism could contribute to the BO 
spin-lattice relaxation in sodium nitrate at room 
temperature. It was shown by Snyder and Hughes (1971) 
that the M-tensor components for this mechanism can be 
expressed as 


es 7 py Sy spy et 
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where Oe is a time dependent field gradient 
component, the bar signifies an ensemble average, and 
ie is the correlation time characterizing the diffusion. 
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we see from equation (7.2) that the M-tensor components 
governing W, are approximately one quarter of those 


2 


governing W Diffusion narrowing of the oN resonance 


Ll’ 
in sodium nitrate occurs near 500 K (Andrew et al, 1962). 
This implies that the correlation time of self- 
diffusion of the Ee ions is ‘olor seconds at this 
temperature. While the correlation time deduced from 
motional narrowing is not necessarily identical to Le 
which occurs in eq. (7.2), it seems clear that eq. (7.4) 


is well satisfied for sodium nitrate at room temperature. 


since the mean deviation, defined in*eq. (Gro eaes 
OC LS*2.0°25 


for the 6 dependence of Wo/W, it follows that the 
"diffusional" contribution is no more than about 0.4% 


t= 1 
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We now compare the experimental M'-values with 

the predictions of the point charge model described in 

Appendix II. Figs. 7.1 to 7.3 show the M'-components 
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Figures aL 


Point charge model estimate of My 3373/1111 
for a "linear" process. The horizontal 


line is the experimentally obtained value. 
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Figure 7.2 


Point charge model estimate of M3333/Mq111 


for a "linear" process. The horizontal. line 


is the experimentally obtained value. 
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Figure 7.3 


Point charge model estimate of My 433/"1111 
for a "linear" process. The horizontal 


line is the experimentally obtained value. 
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Vibration. It can be seen from these diagrams that 
the M' values are determined primarily by 0 rather 
than k. In sodium nitrate the angle 0 is 54.5° at 
room temperature (Cherin et al, 1967). From the 
mMewerOnwCleerace1 Onnstudysorsraul and Pryor (1971) we 
deduce that the parameter k is roughly 1.05. Thus, 
we can see that for M! 


iRS yes 


between theory and experiment is surprisingly good. 


! 
and M3333 the agreement 


aeever unless we assume values of k which differ 
greatly from unity, agreement for M1113 1SeDOOT wen 
Figs. 7.4 to 7.6-~are shown the M'-components for a 
"quadratic" process. Again the agreement is excellent 
for ye eae and Magee: but we now obtain good agreement 
LOy i aac 


sodium nitrate the dominant relaxation mechanism is a 


THUS ,) tte secempting tComconclude that in 


"quadratic" one at room temperature. (This presumably 
would be the harmonic Raman process considered by Van 
Kranendonk (1954), since the anharmonic process (Van 
Kranendonk and Walker, 1967, 1968) is "linear. ) LhasS 
point will be discussed further in section 7.4 when we 
discuss data taken at 77 K. 

The theoretical estimates of Van Kranendonk and 
Walker (1967, 1968) indicate that the anharmonic Raman 
process is roughly two orders of magnitude more 


important than the harmonic Raman process. 
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Figure 7.4 


Point charge model estimate of M1319 "1111 
for a "quadratic" process. The horizontal 


line is the experimentally obtained value. 
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Figure 7.6 


Point charge model estimate of My913/ "1111 
for a "quadratic ™ process. sic horizontal 


line is the experimentally obtained value. 
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However, subsequent more detailed calculations by 
Pietila (1974) have shown that the two processes may 
be equally important for 22N8 in sodium chloride at 
room temperature. Thus, the theoretical situation 


is by no means clear. 


7.3 Comparison with previous room temperature results 


The original enhancement measurements on sodium 
nitrate were performed by Pound (1950). The precision 
of his measurements was sufficient to demonstrate that 
a quadrupole spin-lattice relaxation mechanism was 
dominant in sodium nitrate. Furthermore, he showed 
that W/W, was the order of unity at the crystal orien- 
tations ds= 0%, 

Sarnatskii and Shutilov (1972) have performed 
a double resonance experiment on sodium nitrate by 
using acoustic saturation of a pair of energy levels 
separated by Am = +2. The resonance enhancement was 
measured by a steady-state technique using a Robinson 
spectrometer. They also performed a pulsed NMR 
experiment in order to obtain an estimate of Wi- 
Sarnatskii and Shutilov found that 


= ar 
W/W, = 1.6+0.4 
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W, = 0.066 sec (7.5) 
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atecnesCrystalOnrrentration.le= 50% or 150°%, (The 
authors did not distinguish between these two values; 
néither did they quote the angle $¢.) Sarnatskii and 
Shutilov estimated the paramagnetic impurity concen- 
tration in their sample to be Vl part in Wi on the 
basis of electron paramagnetic resonance measurements. 
By inspection of Fig. 6.5 it can be seen that tolerable 
agreement for W,/W, is obtained if we assume that their 
data were obtained at @ = 150° and > near zero. 
Similarly, their estimation of Wy is in fair agreement 
Pich ourevalucsote (a0 s2r 0084) x 100- sit. Lastly, 
their estimate of W3 is an order of magnitude larger 
than our value (see eq. (6.11)). However, this is 
reasonable since our paramagentic impurity concentra- 
tion is believed to be less than 1 part in Too% 
Reed (1970) measured the orientation dependence 
of W/W, of Ne in sodium nitrate at room temperature. 
It is of interest to compare these results with ours 
Saver t heed ..00/.), eSitice DOCieSstUdITeS were made jon 


the same sample. The ratios of the M-components given 


by Reed (1970) and Hughes and Reed (1971) are 


My 3343/1111 =S 0 so ooje0 6015 
M3333/M@1411 = 0.8222 0.016 eG 
My 433/"1111 =-0.120 + 0.018 
M sO, 142 0,017 
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We can see that the results are in poor agree- 
ment, especially for tg Ava which we found to be 
zero and which should be zero on grounds of crystal 
symmetry. 

In order to explain the difference between the 
two sets of data, we refer first to the 0 dependence 
of Wo/Wy d¥splayedsinetig-s0,o.° sin this figure 
the dotted line is the least squares fit to Reed's 
data. It is obvious that there is quite a large dis- 
erepancy especially near.o = 0° and 180°, Although 
Reed employed a different detection method as 
described in section 9.1, the results presented in 
section 6.3 indicate that the discrepancy cannot be 
attributed to this. In view of the considerable diffi- 
culty we experienced in eliminating systematic errors 
from our measurement system, we are led to consider the 
possibility of some undetected systematic errors in 
Reed's measurements. This suspicion is reinforced by 
the following evidence. Some Wo/Wy values were 
obtained with our system before the elimination of 
the systematic errors caused by the interaction of the 
rf and modulation-related signals. Those obtained near 
® = 0° were in fairly good agreement with those of 
Reed! 

The discrepancy in the Wo/W, measurements is 


displayed in a slightly different form in Fig. 7.7. 
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Figure, /2/ 


Discrepancy between Reed's W/W, data and 
the present data, as a function of the 
difference between the observing and 


saturation frequencies used by Reed. 
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The horizontal axis is the difference between the 
observing and saturation frequencies used by Reed. 
(The negative frequencies in Fig. 7.7 are because 
measurements were made on the Se Ean: eaceleite 
rather than on the "high frequency" satellite. 

From eq. (3.1), it can be seen that near 6 = 90° 
this resonance has a frequency lower than that of 
the center line.) In Fig. 7.7 the data points are 
those of Reed, and the solid line is a linear fit to 
the data. From this figure it seems that the discre- 
pancy is proportional to the frequency difference. 

We shall now describe two possible explanations 
for the discrepancy in the 6 dependence. The first 
explanation presupposes a systematic shift in the 
baseline. We assume that in Reed's spectrometer there 
was a modulation pickup or a radiation problem of the 
type discussed in section 5.4.3. This would cause a 
baseline offset at the output of the lock-in detector. 
We further assume that this baseline was of the form 
shown in Fig. 7.8, where the effect has been exagge- 
rated for demonstration purposes. If account was not 
taken of such a baseline offset, then the measured 
resonance amplitudes near 98 = 0° would be too large 
and the measured value of E would therefore be too 
small. Similarly, near 9 = 90°7 the measured resonance 


amplitudes would be too small, leading to values of E 
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Figure 7.8 


Explanation of the discrepancy in the 


W/W, data in terms of a baseline offset. 


PAIR 


O9l 08 


[ere 3 00=@ 4 


24!]|240S 
Aduenbesy ybipy 


00629 10 
AUG fee 3215 
Aduenbei4 MO7 


euljesog pewnssyy ——— 
QUIOSOG |DOY ceccceee 












a 
; ae 4 
j é , > .. 8 
n : f t . te 
» : . 1 =! 
n Gs 
"“ Mm ' ve) a 
é a 
. ' j “ ” re 


>. 0 


.--yonetpett sty ; 
= > —gilieies. 9 = | 


tt s 7 = o 9 . 9 f j m ; 
* ‘ ; 
“4 . ' . ‘ ¥ 


e 





which were too high. According to this interpretation, 
it is the difference in the baseline offset between one 


peak of the absorption derivative and the other peak 


which causes the error. To explain the linear behaviour 


ee 


Shown in Fig. 7.7, we need a baseline offset of the LOrm 


B= B, + B, (Av)? : na 


where Bo is arbitrary and Av is the difference between 


the observing and saturation frequencies. To explain 


the magnitude of the discrepancy in Wo/Wy the difference 


in the baseline on either side of the resonance corres- 
ponding to Av = 160 kHz, would have been little more 
than the peak-to-peak noise amplitude in an individual 
run. Thus, such a problem could easily have been over- 
looked. 

We shall now consider the second explan- 
ation for the discrepancy. As mentioned in sections 
5.4.3 and 7.1, there was a small frequency dependent 
systematic error in our system, which we believe to be 
caused by the frequency dependence of the rejection of 
the balance network. As can be seen in Fig. 5.4, the 
orthogonal mode of the balance network contains a 
delay line. Thus, since a delay line produces a con- 
stant time delay, the phase angle of the "orthogonal 
mode" will vary linearly with frequency. Thus, if the 


balance network is adjusted to reject the saturation 
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frequency, it will not be perfectly balanced at other 
frequencies. Indeed, ae magnitude of the net coupl- 
ing will be proportional to the difference between 

the saturation frequency and the frequency under con- 
Sideration. Thus, if the rf power source receives 

any radiation from the spectrometer, at the observing 
fPeequency Ani particular pert will ‘be coupled back to 
the spectrometer sample coil via the mechanism just 
described. Since the magnitude of the effect will be 
proportional to Av, this mechanism can account for the 
linear dependence of the discrepancy displayed in 

Fig. 7.7. The magnitude of this problem is likely to 
be more serious in Reed's apparatus than in ours since 
his spectrometer and rf power source were not effec- 
tively shielded. 

Although there undoubtedly was a modulation 
pickup problem in Reed's spectrometer, we feel that 
the second explanation is more appropriate. One reason 
is that the second mechanism will result in a direct 
modulation of the spectrometer observing power, and 
hence produce a systematic error which is independent 
of whether ow, oF 2W detection methods are used. 

Thus, it can explain the agreement between our pre- 
liminary results and those of Reed. Secondly, the 
discrepancy will automatically be zero when Av is zero, 


whereas in the first explanation one has to. make a 
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separate assumption to obtain this feature. 

The discrepancy in the 8 dependence, just dis- 
cussed, is responsible for the difference in the values 
of M 


and M3333/M@1 441° However, the gross 


1313/1111 


discrepancy in the value of M stems from the 


Tea Quy tid 
measurements made by Reed near ¢$ = 90°. His data are 
shown in Fig. 7.9 together with @ dependence obtained 
from our M'-values. It can be seen that the 9 depen- 
dence of W/W, should be almost symmetrical about 

6 = 90°, whereas Reed's data are very asymmetric. 
Moreover, the discrepancy between our values and those 
of Reed does not appear to be proportional to Av,in 
Seouteractutoe tic datassnowniei ny hida av. 7... —Since no other 
data were taken by Reed after these very anomalous 
results, one cannot say whether or not a permanent 
change in the characteristics of the spectrometer 

had occurred after the 8 dependence results. AS was 
made clear in section 5.4.3, apparently insignificant 
changes in our apparatus did produce large systematic 
errors. We can only speculate that some change occurred 


in the equipment between the two sets of data. 


7.4 Comparison with 77 K results 


Goldburg (1959) performed a double resonance 
experiment similar to ours except that he measured the 


enhancement by a pulsed rather than a steady state 
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Figures) .o 


W./Wy as a function of angle 9 for $ = 86.5°. 
The solid line represents the present study, 
and the dotted line is a least squares fit 


of Reed's data to the theoretical form. 
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technique. He found Wo/W, Sr Omg Ube Uo at Ue Ol ator 
sodium nitrate at 77 K. 

Niemela (1967) studied the orientation dependence 
of the eon relaxation in sodium nitrate at 77 K by 
means of the pulsed NMR analog of our transient 
experiment. He fitted the transient decays to a sum 
of exponentials from which he was able to extract 
values of Wis Wo and W3 ac VarlOuSsecrystal orientations. 
At 9 = 0°, Niemela's value of Wo/Wy is approximately 
1.06, in only fair agreement with Goldburg's result. 


The following results were derived by Hughes et al 


(1970) from Niemeld's data: 


TC ht gee See Oe 
M3333/81171 =U OG ter 04 Giz) 
Ma aa3/Miii1 = 0-264 0,02 

and 
e°0°M 434% (1-844 0,05) x1074 sth. (7.9) 


It is known that thermal expansion in sodium 
nitrate occurs mainly along the triad axis (Austin 
and Pierce, 1933). If we assume that the thermal 
expansion coefficients obtained from X-ray measurements 
by Cherin et aly (1967) mapoividown tol77 k,5 then the 


angle 9 in the point charge model will increase by 
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approximately 0.5° on cooling from 295 to 77 K. Thus, 
Erom’ Figs. 7#lland’7s2"we' sée that Mi 373741111 and 
M3333/M14071 should both decrease by approximately 0.03. 
This agrees well with the observed decrease of 

OO 2120 028and 0805 210 504 respectively between 295 and 
77 K. However, we see that My 13/1111 is more than 

a factor of two larger at 77 K than at 295 K. (The 
sign of Niemela's value of My413/ 1111 is not known, 
Since he did not distinguish between d¢ and ¢+ 60°. 

We assume, for lack of evidence to the contrary, that 
Mi333/™@iiii is negative at 77 K.) By comparison of 
Btqsts/.3°anae 736 it is.clear that 'theW77 °K "valuevot 
My443/1111 1S in much better a Graenet with the 
"linear" model than with the "quadratic" model. This 
would suggest that the dominant relaxation mechanism 
is the anharmonic Raman process at this temperature. 

In section 7.2, we saw that our room temperature 
results are consistent with the harmonic Raman process. 
However, according to the theoretical predictions of 
Van Kranendonk and Walker (1967), both processes should 
have the same up? temperature dependence. It thus 
seems unlikely that the dominant relaxation mechanism 
would be different at-295 and 77 K. One is thus led 
to question the validity of the point charge model. 

It seems reasonable to expect this model to be more 


appropriate at high rather than low temperatures, 
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Since, at lower temperatures the characteristic wave- 
length of lattice vibrations is longer. Thus, one 
would expect more distant neighbours to contribute to 
the quadrupole relaxation at the lower temperature. 

On this basis, one is led to conclude that the harmonic 
Raman process may be the dominant one at boths tempera— 


tures for =— 


Na in sodium nitrate. However, this is very 
tentative and more work is required to confirm this 
point. 

BrOMmeces aru acarand © (6a! 2jaswe wind athat eaGe Mao 
has increased by a factor of 26:24 0.8 between 77 and 
299 Ke Accordingito eq (2.40);, fthis «quantity ‘should 
have a 7? dependence and hence change by a factor of 
asia However, since 7/K as probably Jess than half 
the Debye temperature for sodium nitrate, the more rapid 
variation of 270M, 313 May not; be too’ surprising. 


Niemela (1967) found that the magnetic relaxation 


probability is given by 


fey = (a a) hire (Fal) 


at 77 K in his sample. Our value of the magnetic re- 


laxation, 1s given by (fromesect ron, 6.4) 


W, = (2.74 029) a0n en F (Jeti) 


Since magnetic relaxation is caused by paramagnetic im- 


purities is expected to be temperature independent, at 
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least to first order, it would appear tromeqs. (7-10) 
and (7,11) that our small crystal is of comparable 
purity to Niemeld's crystal. Since the paramagnetic 
impurity concentration of both crystals has been 
estimated to be no more than 1 part in Looe the agree- 
ment between the W3 values is reassuring. 

The excellent agreement between the values of 
Mi373/@1111 and M3333/@11071 obtained by us and by 
Niemela is also very encouraging, since the experi- 
mental methods and temperature were very different in 
the two cases. Thus, we feel that the large change in 


. taps 1a: ee 
My333/"1111 Sy Sayles! pe further investigation. 
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CHAPTER VIII 
CONCLUSION 


From the results and discussions presented in 
Chapters VI and VII, it is clear that we have developed 
instrumentation capable of obtaining precise and 
accurate measurements of NMR amplitude. This has 
enabled reliable relaxation data to be obtained on 
the quadrupole split spectrum in sodium nitrate. 

Since the system is under the control of a real-time 

minicomputer, very complex control algorithms can be 

executed. Thus, experiments which would otherwise be 
impractical can be conceived and performed. 

Our steady state method suffers from the dis- 
advantage that it does not give the magnitude of Wy or 
Wo- However, much information on quadrupole relaxation 
can be obtained from the ratio W/W, - We have shown 
that this quantity can be measured using our experi- 
ment and technique to an accuracy of about 41%. 

The transient method has the advantage that it 
measures the relaxation rates directly and is much 
less demanding on the instrumentation than the steady 
state one. It would appear from our preliminary work 
that it is capable of giving values of relaxation 
rates which are accurate to about +1% for the case 


where there are several closely spaced resonances. 
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(In this case, the pulse method suffers from the 
drawback that more than one resonance is excited 
Simultaneously thereby complicating the interpretation 
of the data.) 

We have definitely established that the = gp 
spin-lattice quadrupolar relaxation occurs via a 
multiphonon interaction in sodium nitrate at room 
temperature. By comparing the results with a point- 
charge model, it would appear that the spin-lattice 
relaxation occurs via the harmonic Raman process at 
room temperature. The evidence at 77 K is that the 
harmonic Raman process is not dominant. However, this 
conclusion may be due to the limitations of the model. 
The theoretical and experimental orientation devendence 
eee Wo/W, are in excellent agreement, and no Signor 
the anomalous behaviour reported by Reed (1970) is 
found. Furthermore, we conclude that the anomaly was 
the result of systematic errors of instrumental origin. 

We now present some suggestions for further 
investigation on sodium nitrate. It would be interest- 
ing to investigate the temperature dependence of 
My473/%1111 in view of the large change in this quantity 
far from any known phase transition. Measurements of 
W/W, in the vicinity of the phase transition at 276°C 


should be useful in studying the detailed mechanism of 
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the transition. In particular, since "soft vibrational 
modes" would produce a direct relaxation mechanism, 
we would expect them to have an effect on the value of 


14, in sodium 


W5/W) - A double resonance study of 
nitrate would be useful in providing further evidence 
on the relaxation mechanisms Scecurminguin this crystal. 
Measurements of Wo/W, byallLtmasonic=saturation in 
many cases, gives rather unreliable results. Since 
we have a detailed knowledge of this quantity in sodium 
nitrate, it would be interesting to perform such an 
experiment, using our experimental aiGangencnt co gery 
and determine if there are any fundamental differences 
in the two techniques. It would be useful to extend 
the Wy measurements to other crystal orientations to 
see whether the data can be fitted to the theoretical 
form. This would check whether there are any unsus- 
pected systematic errors in the Wy determination. Also, 
the transient measurement may possibly be extended to 
the measurement of Wo by monitoring the transient 
behaviour of the center line and fitting the decay to 
a multiexponential time dependence (cf. Niemeld, 1967). 
Other sodium salts which would be suitable for 
a double resonance study are sodium nitrite, sodium 
chlorate and sodium bromate, for example. Sodium 
nitrite is particularly appealing since it has several 


phase transitions. Furthermore, since the sodium sites 
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in this crystal do not possess inversion symmetry, 
this material should provide a test of time reversal 
symmetry (Hughes, 1973). 

The work described in this thesis was carried 
out on single crystals. However, there are many sodium 
Salts whichsare difficult, af not ampessible, to obtain 
as a large single crystal. However, if the asymmetry 
parameter is zero and the quadrupole coupling constant is 
sufficiently small (probably $500 kHz), a peak in the 
NMR powder pattern may be visible, corresponding to the 
satellites for which § iS near 90° (Cohen and Reif, 
1957). While this peak. will be weak, it may be possible 
to extract meaningful values of Wo/W, and Wy separately 
for 6 = 90° by the steady-state method and the tran- 
sient method described in this thesis. An obvious 
sample with which to test the method is sodium nitrate, 
since the quadrupole relaxation probabilities are 


already known. 
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Appendix I 


General orientation dependence of Wy and W, 


We give the general orientation dependence of 


the quadrupole relaxation probabilities Wy and Wo, 


the case of a real quadrupole relaxation tensor 


£Or 


(Snyder and Hughes, 1971). 
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Appendix II 


Calculation of Migarg! for sodium nitrate using a 


point charge model. 


We shall describe a simple nearest-neighbour 
point-charge model used by Hughes and Reed (1971) 
to calculate the M-tensor components in sodium nitrate. 
The sodium nuclei are assumed to vibrate in the elec- 
trostatic field of the nearest neighbour oxygen atoms, 
which are treated as stationary point charges. The 
configuration of these six charges is shown in Fig. 
Poet | 

The change in the electric field gradient Noe 


at the nuclear site caused by a displacement dm can 


be expressed by the following power series 





2 
OV a V 
e ap i OB. 
Vag=) a ox; +5 2) Ray axyda, +s. eee, 
ibe Lig Es 


The linear term in eq. (A2.1) will give rise to spin- 
lattice relaxation via the direct process and the 
anharmonic Raman process (Van Kranendonk and Walker, 
1967, 1968), whereas the quadratic term will give rise 
to the harmonic Raman process (Van Kranendonk, 1954). 


The M-tensor components M are governed 
a. 


Ba'B' 
by the product of the fluctuations in Vag and Vargre 


Thus, for a "linear" process. we obtain 
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Figure A2.1 


Schematic representation of the charge 
configuration used in the point charge 


model for the calculation of M ess 
apa's 
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where the bar signifies an average over the nuclear 
vibration. If we assume that there is no preferred 


direction for the nuclear motion, then 


ee 2 
dx dx, = 6., (dr)“/3 . (Oe) 


substituting into eq. (A2.2) we thus obtain 
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(M : (A2.4) 
The model can be generalized to take into account 
sanisotropy in the nuclear motion.i° Let k be the am-— 
plitude of the vibration along the z axis divided by 
the amplitude in the x-y plane. We further assume that 
the vibration in the x-y plane is isotropic and of 


amplitude dr. Then eq. (A2.3) will become 
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For a "quadratic" process, we obtain from 


eq. (A271) 
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for the case of the nuclear site possessing inversion 


Symmetry. For the anisotropic nuclear motion we obtain 


See 4 
4 _ 4 Cadi) 
(dx) =k bE a (A225 7;) 
and 


ee oe 4 
2 ie ape (dr )ay 
After a lengthy calculation, which takes into 
account the configuration of the six charges shown in 


Fig. A2.1, it can be shown (Reed, 1970) that for the 


"linear" process: 
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+ (3/2) (-1+5sin78/4) ]+k*éos76[ (5sin76-1) 7 
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+ k*[3-30cos7e+35cos9}]7) . (A2.15) 


noe Merl te Mees -10 ae 
Mi4137 (34 u sin@cos$cos3$ r /2) [3 (7sin 6-3) 
x (3-30sin79+35sin‘6) -3 (5cos76-4) 
- (3-30sin?6+35sin‘6)+(7sin26-3) (5cos29-4) 


- 3(3-7sin*6/4) [3- (30/4) sin76+ (35/16) sin’6] 
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3 (-14+21sin70/4) [5cos76-4+ (105/16) sin‘6] 


2 


+ k7[(35sin7ecos6-4) (42sin29-20) -16 (7sin26-12) 


x (15sin“6-16+35sin~6cos“0) -60sin~6 (7cos29-1) 


(Isin@6=1)\(7sinse24)-9 (8-7608-26) 


x 


* (15sin76-16435sin~-6cos-6)] 


9 
tien 12) 7eose6= 3 (oSseincecossesisei na niie. © e(aola6) 


In the coordinate system and for this model Mi 193 aS 
zero. 

Since the Mga! are proportional to gos fon 
the "linear" process and to po Bee the "quadratic" 
process, the nearest-neighbour approximation is likely 
to be reasonable for the case of a ine sOdium nitrates 
Also, although the details of the phonon spectrum have 
been ignored, the model should be a reasonable 
approximation at temperatures greater than the Debye 
temperature. The point-charge model is a gross simpli- 
ELCactOneor thes actuallsituation ingarrealscrysctaln 
Moreover, factors such as antishielding (Sternheimer, 
1950) and covalency in the crystal bonding (Yosida and 
Moriya, 1956) have been ignored. This will result in 


the magnitude of an individual M component being com- 


pletely unreliable. However, the ratio of two M 
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components should be fairly accurate, since it is 
expected that several factors which have been ignored 


will affect all the M's equally. 
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Appendix III 


Selection of Oe values and the optimum utilization 


of time in the measurement of NMR enhancement 


In our experiment, the NMR enhancement E is 


measured at various observing powers P and the 


obs’ 
value of W/W, is determined from the intercept on 

the E axis of the straight line which is the best fit 

to the E-Pobs data. The enhancement values are subject 
to random error because of the noise generated in the 
spectrometer. However, for a Robinson spectrometer, 
this noise is independent of rf level (Howling, 1966; 
Smith and Hughes, 1971). At the low spectrometer powers 
used in our work, the NMR saturation is negligible and 
the signal amplitudes are therefore proportional to the 
rf level Eee Since E is only a weak function of 
ADs (it only varies by 5% over the complete range of 


SDE used), it follows that the random error in the i'th 


data point Ey, is inversely proportional to Y(P ope 
Also, the error in E; is inversely proportional to vti, 
where t, is the-time spent) on the i’th datavpoint.. On 


ee and 


this basis, we wish to find the values of (Pore : 


t; which minimize the random error in the intercept 


subject to the constraint 
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where T is the total time available. Since many hours 
are required to obtain a. precision of 0.1% in an en- 
hancement measurement at even a large value of Peas 


it is important to acquire data as efficiently as 


possible, 


If we now call (P and Eis Xs and y; res- 


anes 
pectively, then the data are fitted to 


y=a+t bx (A3e2) 


where x is necessarily positive. The square of the 
StanGabdeecrrol in the intercept a is given by 


(Bevington, 1969) 


oF = x2 /o2/{ (J1/04) (1x4 /o2) - ()x,/o5)7} (A3.3) 


where oO. is the standard error in Ne and the summations 
extend over all n data points. It follows from the 


discussion in the previous paragraph that 
NV Moi  fapienien. (A3.4) 


where c is a constant. If we regard the xX, as arbitrary, 


then the optimum t, must satisfy the relation 
d(o7) = y(307/at Vdter=20 (A325) 
a a al a ; 


subject to the constraint 


Wels Se (A3.6) 
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A solution to eqs. (A3.5) and (A3.6) is 
2 : : 
a0./at. = a constant independent of 4 . (AS) 


Substituting eq. (A3.4) into eq. (A3.3) we find 


30° i ee 

x= Se ee ee (A3.8) 
; S 2 22s , 
j AVE Shore thee Gene) ) 


Since the denominator in eq. (A3.8) is independent of 


pele LOlLiows fromeeq. (42.7) that 


. 3 ‘ 2 Zan : 
ete x)xit,) a constant independent of j. (A3.9) 


since eq. © (A323. %).ts taecubic tequation .in a there are 
at most three independent x values which satisfy eq. 
(A3.9). Expanding eq. (A3.9) for n = 3, we find that 
one of these equations is a linear combination of the 
other two. We therefore conclude that the optimum 
number of data points is two. 

We now consider the specific case where the two 


data points are at 


mss 
(A3.10) 
Xo = ax ) Wiz Ch OS ec <ae- 
and the corresponding measurement times are 
bere (A3.11) 
toa Cl—B LT monte On es oh 
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Substituting into eq. (A3.9) we find that 
ae (A3.12) 


This gives the optimum allocation of measurement time 
for an arbitrary value of a. To find the optimum 
VeluerOL awe Substituteseqs, (A3.10), (A311) and 


(A3.12) into (A3.3) to obtain 


2 
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By differentiating eq. (A3.13) we find that the optimum 

value of a which minimizes ot ous AP eeeromeed et AG ol 2) 

one finds that the corresponding value of 8 is gS 
Thus, we conclude that the optimum method of 

data acquisition consists of taking just two data points 

avitering by a factor of four in Peter and spending 


eight ninths of the time at the lower data point. 























